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ABSTRACT, 


Conformal transformations are given which are suitable for the discussion 
of the space outside a finned rectangle. These include figures such as 
rectangle, cross, three-sided box, form like the letter H, ‘‘ windmill ’’, 
and the transformation could be applied to any rectangle with four fins 
of any lengths at its corners, each of them in either of two directions. 
The transformations convert the outer space into the inside of a rectangle, 
which can then be transformed back into the outside of a circle, with the 
same point at infinity as the original space. 

By a generalization of a special case of this transformation another is. 
given, whereby an irregular star, composed of any number of rays of any 
lengths and at any angles branching out from a centre, can be trans- 
formed into a circle. This includes, as a simple example, a “ dihedral ”’ 
of two equal arms at any angle. 


1. THERE are many uses to which conformal transformations are put 
in physical problems, and the discovery of the appropriate transformations 
is not always easy. The only general systematic method of discovering 
them is the Schwartz Christoffel method, and this is limited to cases. 


* Communicated by, the Author. 
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where only half the plane is of interest, and to cases where this half-plane 
can by symmetry be extended to describe the whole plane. As soon as 
less symmetrical shapes are considered, the discovery of the appropriate 
transformations becomes guess-work, and it may therefore be useful to 
put on record a set of such transformations representing a rather wide 
range of forms. 

The present note started with some work of V. M. Falkner* in connection 
with ‘finned aeroplanes, for which I provided him with certain trans-— 
formations involving elliptic functions. For this work the Schwartzian 
method sufficed, but I have had the curiosity to extend it further into 
a range of shapes for which that method cannot be used. It is hardly 
necessary to say that, though it was aerodynamics that gave rise to these 
problems, they are not to be interpreted as giving directly the air flow 
over an aeroplane of the shape described. 

Elliptic functions will rationalize four square roots, and therefore any 
form involving four right-angled turns can be represented by their means. 
There is, however, one important further limitation on the shapes here 
considered, in that I shall only consider the transformations mapping 
the area outside the form and extending to infinity in all directions. 


Fig. 1. Fig. 2. 


The present work gives the transformations whereby the area outside 
any finned rectangle like that in fig. 1 can be transformed into the inside 
of a rectangle. This rectangle can then in turn be transformed into the 
area outside a circle in such a manner that the point at infinity coincides 
with the original point at infinity. In the course of this work the con- 
sideration of one special case, the cross, suggested another generalization 
by which it is possible to transform into a circle a figure that may be 
described as an irregular star, composed of any number of rays of any 
lengths and at any angles, such as shown in fig. 2. Perhaps the most 
useful example of this would be an * aeroplane with dihedral angle ”. 

2. Since the purpose of this note is chiefly to put the transformations 
on record, there is no need to explain in detail how they were found. It 
need only be said that there are two primitive forms from which they were 


* Aer, Res. Co., Rep. & Mem. No, 2280 (Oct. 1945). 
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generalized which may be called the H and the X, from the shapes of 
those letters. The H was the subject of the previous note, though the 
notation has been improved from that used there. The interesting case 
is where the cross-bar is not exactly half way up it. By bisecting the H 
vertically the Schwartzian method can be used, and the most important 
point that emerges is that if it is to fit together properly, the elliptic 
integrals of the third kind which arise in the integration must vanish. 
This is then seen to be a general principle, because such an integral has 
a logarithmic singularity, whereas the only admissible singularity is a 
simple pole to represent the point at infinity. This condition very much 
simplifies all the work. Other problems involving four right-angled turns 
can also be solved by elliptic integrals, but in general they involve integrals 
of the third kind. 

‘The X, or cross, does not lend itself directly to the Schwartzian method, 
though a symmetrical cross can be discussed in that way by twisting it 
through 45°. However, as will appear, it is very easy to guess the 
transformation directly, and there is therefore no need to explain this 
process. 


—K=<K" 


3. In order to follow the various diagrams it will be well to explain 
their lettering. The four right-angled turns are lettered ABCD taken 
counter-clockwise in the z-plane. Since the outer region in z is to be 
transformed into the inside of a rectangle in ¢, these points go clockwise 
in ¢ as shown in fig. 3. The corners of the rectangle are the points 
4+K-+iK’ in the ¢-plane. Any point on AB is labelled Q, and similarly 
for the other sides by R, Sand'T. In particular these letters will denote 
the turning points of the fins; if there are more than one on a side of the 
rectangle, they are given subscripts, and the subscript 0 will be used in the 
particular case where the mid-point of the side is referred to. One other 
place of importance in the ¢-rectangle is the pole P, which corresponds 
in the z-plane to the point at infinity. er 

4. The transformations all involve expressing z in terms of elliptic 
functions of €. In addition to the ordinary functions, we shall have elliptic 
integrals of the second kind, and the most convenient one to use 1s the one 


B2 


4 Sir Charles Darwin on Some Conformal 


that vanishes at both 0 and K. In the ordinary notation this is written 
with a term in E/K, but in the present work E will never occur by itself, 
so we shall adopt a single symbol and write 


CN 8 Ad Aled Sgt ag eae # id (4.1). 
Then take ; 


¢ 
C)=| dn? w du—Gt. oats Gare eS Ea 


This function, like sn¢ etc., is infinite at the points +7K’, and this is. 
inconvenient. It is quite easy to deal with the infinity of sn @, ete., but 
it is better to eliminate it for Z by considering the function 


Zity SECT ORNs terme 


sné 


instead. The second term destroys the infinities at +iK’, but makes. 
a pole at the origin. It is then comparatively easy to shift this pole 
to any other point P. 

To work out the values in z corresponding to the rectangle in ¢ it is. 
convenient to have a table, for though the values of the functions are all 
given by well known formule, there are a good many such formule, 
and their signs are rather confusing. Table I. gives the expressions used ; 
the modulus is taken as & for sn, en, dn, Z, but as k’, its complement, 
when they are accented, as in sn’, en’, dn’, Z’. The last row gives the 
ordinary Jacobian transformation. 


TABLE I. 
C sn en dn ¢ si ecoe C/sn € 
é 1 tk'sn'n | k'en'n {zie c 
K pee ee Dt pasa a es Ss , 
Q +1 dn’ dn’7 dn’ v Z'(n int st t 
oy 1 idn é icon € lar 
ENGI teeta ksn & ksn€é sn. & Ho) 3K 
] ik’sn'n k’ cen'y {2 
ed PEA vege ety Paleake pas 
Bom dn'n dn'n dn’ 9 4) + sR KR sex} 
aaa 1 * 7dné icon é qa 
ale —— ee ae, age ES, 
sae 3 ksn & ksn& sn & A(¢) 2K 
ah en's 1 Hh =i 2) EN aT } | 
| | cn’ + cn’ + en! Tr sn! 2KK"| | 


The following points may be noted. On the boundary sn ¢ is everywhere 
real, being always greater than 1 and positive on the right-hand side, 
negative on the left. cn is everywhere purely imaginary on the boundary, 
and its real part is everywhere positive inside. dn is real on Q and 
8, imaginary on R and T, while its real part is everywhere positive inside. 


se 
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The function in the last column is imaginary on Q and §, vanishing at 
the mid-points, while on R and T it has constant imaginary parts equal 
and opposite in sign. 

5. The transformation yielding the symmetrical cross is 


eran ¢ nes 
| Uta Gat (5.1) 
This is easily generalized to an irregular cross as 
d 
of as aaa (5.2) 
sn€—f—igenl 


where f and g are real, since the denominator is then real all round the 
contour. For the pole to fall inside the rectangle it is necessary that 


eta eesl AB Ge toieF ls 44) 8 Se wat 85 -(B.8) 
The transformation yielding the H with unequal ends is 
be dn ené : 
des Aiea (5.4) 


A corresponding generalization for H can be made, and combining the 
two formulae the final result is 


Zz Pyare sn f+ 7th) 


i AS sn €—f—igen¢ 


where all the parameters J, k, f, g, h, m are real. r 

The first two terms are a generalization of those given in the last column 
of Table I. By adding and subtracting dn cn ¢/ sn in (5.5) it is easy to 
see that the second term still destroys the infinity of Z at Ry and Ty, 
and now yields a single pole determined by the denominator of the second 
term. The general formulae corresponding to the last column of Table I. 
are not written out here, since they are rather complicated, though they 
may be worked out easily by the use of Table I. 

The equation (5.5) has six parameters /, k, f, g, h and m, and by suitable 
choices these can be adapted to give the six dimensions of the finned 
rectangle, that is, its two sides and the lengths (and directions) of the 
four fins. 

The turning points of the fins are given by equating to zero 


1 dz 


+m, ay et eros) 


m—G-+ {—(1—f2)-+sn? (1 —h2f2+ bg?) -+-snk'*(—f+-gh) 


~ ag : 
+éenl(—h-+fy)+é sn onf{—g(1-+k"2) +fhk?]}/ (snl —f—ig en 6), 
(5.6) 
which is a quartic in sn. 
Writing p= 1/(1—f?—g") the pole P is at Cy, where 
+igp p—ug 
on = 1-9 5 Chtye= 19 4 5 5 4 * 5 (5.7) 
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6. We now describe some of the shapes given by specializing the para- 
meters. The diagrams are only sketches, in that when it is asserted 
that one parameter is changing no attention has been given to keeping 
the others strictly constant. To draw them accurately would involve 
much work and give little further information. 

Take first the shapes that are symmetrical about two central lines at 
right angles. For these f, g and h vanish. 

The first to be considered is 


2=Z(C)+ 
Fig. 4 


R 
BG a | 
DIA | 
Raion a Thy ap Th 
As k increases, this gives the sequence of fig. 4. The turning points are 
the solutions of 


dnl ené 
snl 


(6.1) 


or transforming to the lines R or T by Table I., 


dnt $= G56.) ee en eee ane 
With this value of €, 
in 
pe yf at fats). talbio aeeSy aeia Ged 
Fig. 5 


Qa Q2 

ek 

sipramary 
Q, S3 Q, 
at the ends of the fins. When k is not large € is somewhere near 1K, but 
tables of elliptic functions must be used to evaluate it accurately. The: 
value of k must be very close indeed to 1 to get an H with long arms 


and a short cross-bar. 
Next take 


=. 

ll : ces | 
r, [>) i Las i To Hs Th D A Ss 
2 


Pa Hat Guanes eet gene (6.5): 


sn¢ 


and consider how it behaves as m changes. The effect of m is to thicken. 
up the rectangle as in fig. 5. It must be positive, for a negative m would. 
cause AB to overlap CD the wrong way round. 
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_ The rectangle has no fins when m=G—k”? after which the fins switch 
through a right angle. At m=7/2KK’, BC coincides with AD, and 
the sequence stops. The second half of the sequence can be directly 
derived from the first by the Jacobian transformation. 

For the special case k=1/./2, m=G—}=n/4K2 the shape is a square, 
and the Jacobian transformation shows its tetragonal symmetry. 

7. Ifnow aterm in h is included, the two lines of symmetry are destroyed, 
but the centre is still a centre of symmetry. 

First take 
dnfen¢€ ithdné 74 
Fear a ce Se OL i an WED 


As hf increases from zero the fins alter as shown in fig. 6. When 
h=(G—k”)/kk’ one pair of fins disappears. As h grows to infinity the 
side-arms close down, leading to a pure cross when h becomes infinite. 
If a term in m is added, this merely thickens up the rectangle. 


z=Z(t)+ 


Fig. 6. 
Q Q 
R; ar Row. Ry al R, R B R B : 
ere 
i 
A D 
A a 

TA | oe T D 

Ss S 


Attention may be drawn to the “ windmill ’’, fig. 7: 
dnfenZ 7 GTN a 
ie Ode nian fe me) 
which is tetragonally symmetrical. For those familiar with the 
peculiarities of the hydrodynamic fluid it is hardly necessary to say that 
no matter how strongly the wind blows, or in what direction, the windmill 
experiences no resultant couple at all ! 
8. Next consider cases where P is no longer at the origin. Take 
dn@enZ 
snf(—f °° 


w=Z(C)+ 


z=2Z(C)+ (8.1) 
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For growing f the family of shapes is then as in fig. 8. The transitional, 
case with only two fins occurs for 


$= (GEIR, oe a eae 


If a term in m is added, the rectangle is merely thickened as before. 
In both these cases a quadratic in sn £ determines the values of € or 7 at 
the turning points, the two pairs of fins being associated with the two 
signs of 7 en ¢. 


Fig. 8. 
a Qa Qa 
R—c 1g Sie ALE ok c R . B 
R 
D T : . , 
emcee ft Mtr) ” “la, Q, 

If further a term in / is added, the matter becomes more complicated 
because the equation determining the turning points is an irresoluble 
quartic, so that numerical solutions would have to be done. Fig. 9 
sketches the shapes as / increases from zero to large values. 

9. The case where f vanishes, but not gy, can be reduced by the Jacobian 
transformation to the case of gy vanishing and not f. When none of the 
constants vanish the matter becomes very complicated. If numerical 
values are assigned to all the constants, it is comparatively easy to find 
the shape. In doing so it is easiest not to use (5.6) in order to find the 
turning points, but instead merely to estimate the maxima from (5.5) 


by taking out. a few values in their neighbourhoods. Thus fig. 1 is a 
somewhat exaggerated sketch given by the values 


Fig. 9. 


l=1, k=0-5, f=0-4, g=0-3, h=0-1, m=0-2. 
The corners are at +-0-337-+70-501 and the fins at 

Q=0-337—1 0-565 R=—0:339+70-501 

Q.=0:337-+7 0-866 T=—0:386—70-501 


The converse problem of finding the parameters for a given shape 
is much more troublesome. I explored the method of solution by taking 
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the comparatively simple problem of finding the parameters for a -cross 
with arms of lengths 1, 2, 3, 4. This has the transformation 

x th dnl 

~ snl—f—ig en" 
Starting values were guessed for k, f, g by considering the values for the 
half-symmetrical cases, where one pair of opposite arms are equal while 
the other pair have the required ratio, and with these the arms were 
calculated. It was also possible to set down an equation for the variation 
of the parameters, so that the effect of Sf, etc., on all four arms could be 
seen; most of the quantities for these variations occur already in the 
main calculation. Three stages led to the values 


h=2-835, k=0-853, f=0-489, g=0-294, 


which are correct to about one figure in the last place. The ends of the 
arms are then at 4, 37, —2, —i. 

A similar method should work in the general case. The corners of the 
rectangle immediately give two relations between /, m and k. For the 
rest the method of trial must be used, but there is one control on the 
starting values of f, ete., which would help considerably. The value 
of dz/d¢ takes a fairly simple form at each of the corners, by the 
consideration that in differentiating (5.5) dn vanishes at the corners. 
‘Then according as dz/d¢ changes sign from one corner to the next or not, 
there will be one fin on that side or else either two or none. In this way 
some guidance can be given as to appropriate starting values for kh, f, g, h, 
and the variation of parameters can then be used to get improvement. 
It may be noted that variation of & is a troublesome process, even if 
sn ¢ instead of ¢ is used as the indeperident variable, so that for this one 
parameter it might be simplest in practice to omit this variation, and 
merely to find the shape twice over for two different close values of kh. 
However, there is no doubt that the solution would be arduous, and it 
seemed hardly worth while to make up a purely artificial problem in 
order to solve it. 

10. It will often be convenient to apply a second transformation, 
converting the ¢-rectangle into the exterior of a circle in such a manner 
that the point at infinity becomes the old point at infinity of the z-plane. 

A simple case will show the method. The symmetrical cross is given by 


—— gs CG san AB, GAC 48 nota (10.1) 
sn C 
Write 
1 
ge ae. eee ey ee (1 (5) 


Then the point at infinity in z) is the same as in z. On the sides of the 
rectangle sn £isreal and cn imaginary, so that | 2 21, and so the rectangle 
becomes a circle of unit radius in the z,-plane. 
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In the general case the transformation is — 


ae 1 (1—g?)—(f—tgp)sn 6+ ( (p+) ores ; (10.3), 
eae begs) sn6—f—ig en¢ 


where p= 1/(1—f?—g"). man 
11. If it is desired to expose the finned rectangle to a “‘ wind ”, it is. 
necessary to form a function that has a pole at P, and that is real on the 


rectangle. Such a function is 


ay AsnC6+i7BenZ (11.1). 
Ae anlesiea ene EF atta Hi : 


in which A and B are real. From this the components of the wind at 
infinity are given by the real and imaginary parts of : 


dw\ Asnl,+7B enf, (11.2) 
dz),  dn&,(enZ,+igsnf)>+th)? ~ ~ ° Be: 


where ¢, is given by (5.7), dn ¢) having its real part positive. 
An irrotational circulation round the rectangle is given by 


w=i0ln2, 


where %, is given by (10.3), since Inzy is pure imaginary on the rectangle. 

12. If it is desired to transform the finned rectangle into a circle, the 
intermediate use of the rectangle is in general indispensable, but there 
is an exception in the case of the cross, since the transformation then 
involves no elliptic integrals of the second kind. For example, between 
(10.1) and (10.2) € may be pei Taking the special case of the. 
equal armed cross with k=1//2, the eliminant is 


= 4222-8) So 
This may be compared with z=4(z)+2)+), which transforms the circle 
into a line bounded at both ends. This suggests immediately a 
generalization to a star of six, eight, etc., equal rays, but a very much. 
wider generalization than this can be made. 
Take 
z=lzy *(z—e%)s(z)— 82)... (2g—etn)’m. . 2. (12.2). 
where 


ALAS nee ee 2 en a cen errr 

To study the transformation put ze, and consider the part where. 
%<@<a,. For this part, using (12.3), 

w=4l exp {i} 2a, im(A,—1)} sin” 4(0—a,) sin” 4(ay—6) . . . sin’» 4(«,—8),. 

(12.4) 

so that z traces a ray in the direction $21 a,A,+7(A,—1), vanishing at 


8=a, and «,, and so having a maximum in between. When a,<@<a, 
there will be another ray in the direction AD a,Astm(Ag+A,—1). Thus. 


oe 
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. there are rays of various lengths at successive angles 7A,, 7A,,.... - TA,, to: 
one another. By suitable choice of J, »—1 of the \’s and the » «’s any 
arbitrary irregular star of the type shown in fig. 2 can be transformed 
into a circle. 

In finding this transformation the \’s can be given at once, but the 
x8 would have to be found by trial. It may be noted that if two adjacent 
% s are taken equal, one of the rays disappears, and so.it may be concluded. 
that if two «’s are close together, the ray between them. will be short.. 
We will be content with two examples. ; 

Take 


Cereal (ens yr eee ay ine LALO) 


This represents a dihedral with equal arms at angle 27f apart, as in fig. 10. 


The maxima occur at zg=exp(+76,), where cos @,=1—2f8. The rays are: 
of length 416°(1—B)!-8. 


Fig. 10. Fig. 11. 
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As a second example I arbitrarily chose the problem of an arrow-head,. 
in which the barbs, at 30° to the shaft, were only half as long, as in fig. 11. 
This is represented by the transformation 
25 1(z)—e!*)¥8(zy 4-1) °F (2956), ie (12.6) 


and withou much trouble it was found that «71° 46’. The shaft is 
then of length 3-348 and the barbs 1-673. 
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Il. The Propagation of Sound Waves in an Open-ended Channel. 


By W. CHESTER, 
Department of Mathematics, Leeds University *. 


[Received August 22, 1949.] 


SUMMARY. 


An exact solution is obtained for the problem of a harmonic wave of 
sound approaching the open end of a two-dimensional channel, and. 
suffering reflection and transmission. Within a certain wavelength. 


La a ak mea ee Mg ens ee oA 


* Communicated by Professor S. Goldstein, F.R.S. 
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range, the wave returning along the channel becomes sensibly plane at 7 
large distances from the mouth, the reflection coefficient being 
-exp [—27b/A], where 2) is the width of the channel and Ais the wavelength. 
Explicit formule are also obtained for the change of phase in the returning 
wave, and for the transmitted energy. 

It is also shown that the solution can be used to deduce the energy 
absorbed from a plane wave external to the channel and approaching 
from an arbitrary direction. 


INTRODUCTION. 


WHEN a plane harmonic wave of sound is incident upon a two-dimensional 

channel, formed by two semi-infinite parallel planes, diffraction effects 
will occur at the open end. Ultimately, a periodically steady state is 
reached, and the object of the present paper is to determine the asymptotic 
form of the resulting periodically steady potential field, assuming that 
the wavelength lies within a certain range so that only the dominant 
mode is propagatéd inside the channel. There is no theoretical difficulty 
in extending the analysis to include a wider wavelength range, but, 
from a physical view-point, the present range is considered sufficient. 

Problems of this type were first considered by Rayleigh using 
approximate methods. For the particular case of an incident wave 
originating inside a circular pipe, Rayleigh (1940) gave an estimate of 
the end correction when the wavelength was large, and also deduced 
that it increased indefinitely with the wavelength for a two-dimensional 
channel (Rayleigh 1904). (The end correction is the correction to the 
length of the pipe in order that the open end may be treated as a loop). 

More recently, Levine and Schwinger (1948) gave an exact solution 
for the circular pipe, valid for a wave-length range in which only dominant 
mode propagation occurs. The present investigation gives the analogous 
solution for the two-dimensional channel, and is inspired by the analysis 
of the above authors. 

A detailed investigation is given for an incident wave approaching the 
open end from inside the channel. The whole of the asymptotic 
behaviour of the wave, both inside and outside the channel, can be 
‘described in terms of the discontinuity in potential across the channel 
walls, and the problem consists, essentially, of deriving an integral 
equation for this quantity which is susceptible to rigorous, analytical 
solution. Rayleigh’s conclusion is verified and, in general, the final. 
results assume a simple form. The end correction and transmitted energy 
are exhibited graphically in figures 5 and 6. 

It is also shown that, having obtained the solution for the above 
problem, the energy absorbed by the channel from a plane wave 
approaching the open end from an arbitrary direction can then be 
‘deduced immediately with the aid of a simple reciprocity relation which 
correlates the radiating and absorbing properties of the channel. 
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ASYMPTOTIC BEHAVIOUR. 


Since the analysis is eventually confined to the behaviour at large 
distances from the mouth of the channel, we begin with some general 
remarks concerning the asymptotic form assumed by the potential both 
inside and outside the channel. 

The origin of coordinates is taken at the mouth of the channel with the 
negative v-axis along the centre-line. The potential of the incident wave 
is denoted by A exp [ik(«—at)] or exp [—ik{r cos (9—6)) — at} | according 
as it originates inside the channel or por. from a direction 5, 
and the final potentials resulting from these two incident waves are 
denoted by ¢,, 4, respectively. 

Inside the channel, the most general solution of the wave equation 
(V°+-k?)4=0, expanded as a Fourier series in y and satisfying the 
boundary condition of zero normal derivative at y—-Ld, is of the form 


f=A exp [tkv]+B exp [—tka]+ z AA, exp [a(q?—k?)?] 


+B,, exp [—x(q?—k?)*] cos gy+ si tA exp [a(s?—k?)4] 
+B, exp [—x(s*—k?)*]) sin sy, or) ae Bae tel (256 Ea uest( 1.) 


where 


pgitGe-al { nw4 4 


(The factor exp [zkat] here, and in the subsequent analysis, is suppressed). 

The potential ¢, is symmetrical about the centre line of the channel, 
so that, in this case, A,—B,=—0. It is now clear from equation (1), 
that the condition for dominant mode propagation is bk<7. All the 
coefficients B,, must then be zero if ¢, is to remain finite for large negative 
xz, and so we deduce that, as «> — oo inside the channel, 


@ eA OxD [thal exp oaien|, 5... 1 (2). 

Similar remarks apply to ¢, except that A’, is no longer zero. This 
reduces the range of k to kb<n/2; then 

d,~C exp [—ika], Be es SS ee eID 


since there is now no incident wave inside the channel. 
In practice, the above restrictions will usually be satisfied. The 
condition bk <7/2, for example, is equivalent to b/A <4, and for atypical 
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sound frequency of 270 oscillations per second, the wavelength is four 
feet. 

Outside the channel, ¢, will behave, asymptotically, like a source 
-situated at the origin whose amplitude varies with direction. Thus 


exp [tkr iy 


re 


ga~F(9) (4) 


The behaviour of 4, is, however, somewhat more complex. If 
= ~<Oy)<7, it must take account of the incident wave which persists 


in the region —(7—6)) <0, <7, and a wave reflected from the upper wall 
in the region (7—6))<0,<7, where 6, is the polar coordinate referred 
to the upper lip of the channel (fig. 1). It follows that the correction term 
.due to diffraction at the mouth of the channel cannot tend to zero at 
infinity along 6,=-+(7—®,) since the continuity of the potential must 
be maintained in the transition from one region to another. 

To overcome this difficulty, appeal is made to the known solution for 
the diffraction of a harmonic wave by a semi-infinite wall. This problem 
was first solved by Sommerfeld (1894) and later discussed by Carslaw_ 
(1898). In terms of coordinates (r,,@,) about an origin at the edge of 
the wall (fig. 1), a wave exp [—vkr, cos (@;—6))] is diffracted to give 
a potential, 


- 


k\t : ‘ (2r,)# cos 3(0,—0) 
(2) exp [—iz/4] {exp [—v7kr, cos (,—6)1 | exp [tkaA?] dA 


—o 


+exp [—vkr, cos (0,;+65)] [ exp [tkA?] an} «ye ate 
“| (2r,)F cos HA: +%) 
It is easily verified that the normal derivative of this expression is 
zero along the wall (8,—-+7), and since 


(2r)teosk6 am\* : exp [ikr (1 +cos6@)] 7 

ie 4 exp [tka] dA ~ (7) exp [0m/4] + iMphirt cos dd = 0<0<7—e, | 
ties exp [tkr,(1-+-cos 6)] 

Bikirt cosyg 7 16 S07 ‘ 


it follows that the above potential includes the plane wave contributions 
in the relevant intervals, together with a term which behaves like 
exp [vkr]/r? times a function of 6,, except along the lines 6,=+(7—6 aye 
‘This last peculiarity of the solution is essential if the potential isto remain 
continuous in the transition from one region to another. 

In the present problem there will be further diffraction effects in 
addition to those given by (5), due to the presence of the second 
- diffracting wall of the channel. Their contribution to the potential must, 
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however, take the form 9(0,) exp [ikr,]/r! for large r, and the asymptotic 
behaviour of ¢, is, accordingly, 


k\t 
p,~ (=) exp [—tr/4—ckb sin ayHexp [—tkr, cos (0,—6,)] 


(2r,)® cos'}(0; 80) ‘ 
«| exp [2kA*] dA+-exp [—ikr, cos (8,+,) 
ay ikr, a 
x| . exp [1k22] da + g(0,) “74 eae AP eT) 
; (217) cos $(, +9») fe) 
outside the channel. 
The factor exp[—ikb sin 6)] is due to the change in origin. An 


incident wave exp [—ikr cos (6—6,)] becomes 
exp [—ikr, cos (8,—6@))—7zkb sin 45] 


when referred to coordinates (7,,6,) at the upper lip of the channel. 


Strictly, equation (7) has been derived for 7/2<@)<z7 since, for 
0 <0, <z/2, the boundary for the incident wave term in fp isa line parallel 
to 6,=—(7—6,) and passing through the lower lip of the channel. It 
can, however, be shown that the above expression does represent the 
asymptotic behaviour of ¢, for 0<0)<z. 


ReEcrprRociry RELATION. 

In this section we show that the radiation characteristics of the channel 
can be used to determine the energy absorbed from a wave originating 
outside the channel. 

To obtain this reciprocity relation, we apply Green’s Theorem, 


[[ a bo—bV)AS=| {ba FP — obo Fehds (8) 


to the two potentials ¢,, 4,, the domain of integration being that bounded 
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by the curves s,, 85,8, (fig. 2),and ” being the outward normal. The portion. 
of the bounding curve denoted by s, is taken along the upper and lower 
sides of each wall of the channel, while s, is a large circle, centre A, the 
origin of the coordinates (r,,4,). Both s, and s; are considered to be at a 
sufficiently large distance from the origin to enable the asymptotic 
relations for oe ¢, to be used. 

Since both ¢, and ¢, satisfy the equation 


(V?+k?)¢6=0 
in the region considered, and since, further, they have zero normal 
derivatives at the walls, equation (8) reduces to 


[ {4.3 $y Fh as +) tea a) — $y Fh dso. . aie 


Substituting from equations (2) and (3), the first integra] becomes 
. : 0 

~ | [ia exp [tkx]+B exp [—tke]) a (C / exp [—tka])—C exp [—tka] 
J-b 


x 2 (A exp [tkx]+ B exp [— ite) | Bees ten: pee ce Ce 
C= — 
In evaluating the second integral in equation (9) we note, from equation 
(7), that the behaviour of ¢, on s, is equivalent to that of the terms 
exp [—tkr, cos (8,—6))—7ikb sin 6)] and exp [—zkr, cos (8,+0))—ikb sin@y] 
in the intervals —(7—@))+«<0,<a and (7—6))+« <0, <7 respectively 
together with a term which behaves like A(6,) exp [ikr,]/r?, except in 
the intervals |6,+(7—@))|<e. It can be shown that the contribution 
to the integral along s; from the two small intervals |@,+(7—0)) |<e 
is O(c) in the limit 7,— ©, while the contribution from the term 
h(0,) exp [tkr,]/7? is . 
LF, exp [ikr,;+z7kb sin 6,] 0 spe (M0, exp 7) 


it 
rt or rt 


= exp [tkr,] 0 exp [¢kr,+7kb sin 6,] 
h(6,) ay Ge mala ee eed d@,=0;. 9.5 (A Th 
integrated over any range of 6,, since the integrand is identically zero. 

Again, a term exp [—vkr, cos 6,] (say) in ¢, will give a contribution — 
{ne Gs [ikr,+-tkb sin 6, | 2 


Fe >— (exp [—?kr, cos 6,])—exp [tkr, cos 4] 


ed exp [¢kr,-+-ikb sin 6,] 
x a (Fe) spl ban), 1 40, 


rt 


= —ikr}| f (9;) exp [tkb sin 6, --ikr(1—cos 6,)](1+-cos 6,) dO, +0(r-?) 


ot se tnin Paap ee : 0 . : 
2 ik| f(9;) exp [tkb sin 6,] cos 0,/2 SE: J eotfamtyoXP [ikr?] an| 


xd9,+0(r-4), . . (12) 
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A simple integration by parts now shows that this expression tends 
to zero as 7; 0 unless the range of integration includes the point 
Hie al 

Thus, for 0<6)<z, the only contribution to the integral along 8s, 
in the limit as 7; > 0, occurs from the term 


exp [—tkr, cos (0; —0))—tkb sin 65], 
m a small interval surrounding 6,;=0,. The term 
exp [—?kr, cos (8,+0))—ikb sin 45] 


has no contribution, in the limit, since it occurs only for (7—6)) <0, <a 
which exludes the point 6,=—4). 
It follows that 


lim {te Ob — dy *) ds 


==— lim nmi 10, ) exp [2hb(sin 6, —sin 6))+2kr, {1 —cos (8;—6y)'] 
Ty 0 4,—46 
x (1+ cos (0,;—6,)) dé, 


3 
=— lim 2tkr? f(0 0)| exp [27kr, sin? 6/2] cos 6/2 dé 
= 


Ty > 2 


= —4ikf(,) { exp [2ike2] dz 


z 
=—sik (F) Oh on ae) . Ue tetas ait tar crcen cis: Nc (18) 
Finally, equating (10) and (13), we have 
o-(5) ze Ls Sg ana ae (14) 


The limiting case 6)=7 is anomalous. There is now no reflected wave, 
but the term which takes account of this wave in equation (5) does not 
tend to zero for 6j=7 and the whole expression for 4, is asymptotic to 
2 exp [tkr, cos 6,] for 0<@)<z. The asymptotic behaviour of ¢, for, 
an, incident wave, exp [ikr, cos @,] is, In fact, 


oN A ‘ ; (27)! sin 30 : exp [tkr,| 
t~(2) exp [—in/4-+ikr, cos 8,])" : exp [ikA*] dA+-9(6,) ee 
(15) 


which, for large r, behaves, as it should, like exp [¢kr, cos ,] for 0 <6, <7. 
The analogous relation to (14) is then 


; Fexp a siea . 5 


Equations (14) and (16) enable us to deduce the energy absorbed from 


SER. 7, VOL. 41, NO. 312.—JAN. 1950. C 


18 W. Chester on the Propagation of 


a wave approaching from any direction 69, in terms of the asymptotic 
behaviour of a wave originating inside the channel. 


EXPRESSIONS FOR THE REFLECTED AND TRANSMITTED WAVES. 

The analysis is now confined to the asymptotic behaviour of the wave 
originating inside the channel which, as we shall see, can be described 
wholly in terms of the discontinuity in ¢, across the walls of the channel. 

To prove this, we apply Green’s Theorem to the potential ¢,, and a 
plane wave exp [—‘kr cos (9@—6,)]. The domain of integration is that 
shown in fig. 2 except that the centre of the large circle is now taken to 
be 0, the origin of the (r, @) coordinates. 

The surface integral again vanishes and, using the asymptotic expressions. 
for ¢, and the fact that its normal derivative vanishes along the walls, 


of the channel, we find that 
{ E@ awe (exp [—7kr cos (@—6,)]) —exp [—ikr cos (@—6))] 


0 exp [tkr] 


9 7] : Seren 
= be [baly=o ay exp [—tka cos 6)—iky sin 6],,-, da 
0 
ca Ip [Puly--o gre [—ikax cos 6)—tky sin 69],-—» dx 
-- | {(a exp [¢kx]+B exp [—ika]) & exp [—7kx cos 6)—iky sin Oy] 
Sy ba 
— exp [—tkx cos 6)—7ky sin 69] = (A exp [¢ka] 
+B exp [—ite)} ds=0, 0 eee dy el Degree: ee: — ai eda ea 
2=—L 
where 
[Pal=¢a)y10—Pa)y—o and [Paly-2= —[¢aly--p= H(z) (say). * (18) 


The first integral, by a similar argument to that used in the derivation 
of (13) is, in the limit as 7 + oo, 

~— 2hertikif(0,) exp [iz)/4]. ook a a 
Substituting in equation (17), and simplifying the remaining terms,. 
we get 

. 1 . . * . 0 
exp [8i7/4](2hk7)!f(@)= lim E sin @ sin (kb sin 6) [ H(a) 
Lo / —jj. 

isin (kb sin @) 


x exp [—zka cos 0] dx+ ny) 
sin 


{A 1-+-cos 6) exp [—7kL(1—cos 6)] 


~B(1=c0s 6) exp [1KL(1-+e08 6)]} ] Se ates is Wie lS cm (OU 
=F(k cos @) (say). 
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Since the asymptotic behaviour of H(x) will be the same as that of 
¢, inside the channel, reference to equation (2) shows that the oscillating 
terms on the right-hand side of equation (20) are cancelled by equal and 
opposite contributions from the integral, and the limit is uniquely defined.. 

It is convenient, at this point, to assume that k has a small positive 
imaginary part which will eventually be reduced to zero at an appropriate 
point in the analysis. Considered as a function of the complex variable 
¢, F(¢) then becomes, for Imf>Imk, 


B(¢)= (4? — 0?) * sin {b(h2 — C?)*} H(¢), (22) 
where 
-0 
H(¢)= | # BIE OX Pel S10]. 25% ay) od > ies (28) 


and is the Fourier transform of a function equal to H(x) when 2 is negative 
and zero for positive «. For economy of notation the same symbol is. 
used for a function and its transform, the distinction being made by the 
argument of the function. 

We may, therefore, regard (20) as the analytic continuation of F(¢) 
in the interval —Imk<ImZ<Imbk, so that, for ¢=k cos 6 


F(& cos @)=exp [3i7/4](2k7)*f(@)=k sin @ sin (kb sin 6)H(k cos @). (24) 


The asymptotic form of H(z) given by equation (2), shows, with the 
help of (23), that H(¢) has simple poles at C=+k. Further, from (20), (22) 


exp [3i7/4](2k)*f(0)=2i7kbA = —2bk lim (€—k)H(¢)=—2bk Res H(6), 


Ck G— 
exp [3i7/4](2kr)*#f(7)= —2tkbB=2bk lim (€+k)H(¢)=26k ie H(¢), | 
6K 2 6h 
(25) 


and the reflection coefficient, R, becomes 


Res H(é) 
B t=—k : 
(pe ee ee a (26) 
A Res H(Z) 
To describe the transmitted wave, we define the function 
27 Gg 2 a 
G6) 7 £(8) eve ee. (27) 


~ 26, A P—[ BP)’ 


where | f(@) |? measures the energy per unit angle radiated in a direction 
6, and 2b(| A |?—| B |?), the total average energy transmitted from the mouth 
of the channel. The function G(@) is therefore the ratio of the time 
average energy radiated in a direction 6 compared with the average energy 
which would be radiated by a simple source at the origin of the same 
total strength. 

CZ 
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Equations (24), (25) and (26) also show that 


2bk| f (8) | 
G())= — a i 2h a eee re 
= TO FU-TRP) 
1 {ksin 6 sin (kb sin 6)\?| H(k cos @) |? (29) 
~ 2bk ~~‘ | Res H(é) ?—| Res H(¢) k 
t=k t=—k 
which depends only on the function H(¢). 
Finally, in terms of G(#), the amplitude of the absorbed wave, given 
by (14), satisfies the relation 


| CP=G(O/G(0).in oe Sa ee 


The object of the rest of the paper is to obtain explicit expressions for 
the reflection coefficient R, and the function G(@), which reduces the 
analysis, effectively, to the derivation of H(¢), the transform of H(z), 
since the relevant quantities depend on this function. 

We show first that H(x) satisfies an integral equation similar to the 
Wiener-Hopf type (Titchmarsh 1937). The usual method of solution of 
this equation is well adapted to the present problem since it is formulated 
in terms of the transform of H(x), with which we are directly concerned. 


THE INTEGRAL EQUATION FOR H(z). 


The boundary conditions at our disposal in deriving H(x) are those 
of zero normal velocity at the walls of the channel, together with the 
known asymptotic forms of the potential at infinity. Such problems can 
be solved by expressing the potentia] in terms of an appropriate Green’s | 
Function. It is known from Green’s theorem that if ¢ is a solution of 
the equation (V*+k?)46=0, continuous along with its derivatives inside 
a closed curve C, and if J” is another solution, continuous except near the 
point P where [~ —log p, p being the distance from P, then 


1 oovemall 
HP)= 5 (PE —8S) as Pei hse Te 


nm being the outward drawn normal to C. 
If, in vector notation, we denote by r, r’, the position of a general 


point in space, and the particular point P, respectively, then such a 
function satisfying the above conditions is 


(kl r—r’ |)=P(w—2’, y—y')= HP r—r' |), Tee 3a} 


where H{(z) is the Bessel function of thethird kind as defined by Nielson 
(Watson 1944). We note that the asymptotic form of H(z) is, for large z, 


He)~ (=) exp [ile—wi4)] SMe MR ety ier ard is) 


We now take the contour of integration in (31) to be the curve 81; 83, 85 
shown in fig. 2. By virtue of the known asymptotic behaviour of d, 
a 
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and J’, the contributions from s, and s, vanish in the limit, and since 
0¢,,/0n=0 along the walls of the channel, we obtain, finally, the relation 


: 7] 
Oe ire ‘ H(x’) Bal ru r—r'|)y-,—-2(r—r’ De | da’, . (84) 


and the boundary condition that the normal derivative of ¢, should 
vanish at the walls of the channel gives 


g , oe , , , ~ 
tae H(a ou: lr | oma fe yr —L (| Lap Ds | . dx =U O (35) 


ifex <0: 
Equation (35) is the integral equation for H(a). For convenience it 
is written 


0 
| H(x')K(a%—2’) das =0:, 4g <All). 4 ‘ c : (36) 
where : 
K ‘ = : ge It v ( 4 rg 5 { f ST 
ae 27 Oyoy’ (@ YY yy —L(@, YY es |e ee OL) 


In deriving the above equation, no restriction to dominant mode 
propagation was made, and the result holds for an incident wave of 
arbitrary frequency. It is in determining the precise nature of the 
solution that the asymptotic behaviour of H(a) isrequired. In the present 
case the solution is subject to the condition that, for large x, 

H(x)~A exp [7kx]+ B exp [—tkz]. Sel ger gibi Atak 
which is sufficient to determine the function uniquely. 

Analogous results were obtained by Levine and Schwinger (1948) 
for the propagation of sound waves in a circular pipe. It was shown 
by these authors that integral equations of the above type could be 
_ solved by the use of Fourier transforms and their analysis is here adapted 
to the two-dimensional channel. 


SOLUTION OF THE INTEGRAL EQUATION. 
In solving (36), the equation is first written in the more general form 
0 
| H(x’)K(a—2’) dx’=E(x), x>0 
“ L2G = ate 0: 


Pe teat ose COD) 


so that, multiplying each side by exp [—7a], and integrating over all 2, 
we obtain 


( H(x’) da’ |e K(«%—x’) exp [—ila] dx= fe E(x) exp [—ila] da, 
or | aoe 3 
f H(x’) exp [—ila"] dx’ {9 K(z) exp [—2¢z] dz= th K(x) exp[—ila] da, 
fe | 
BOK Ow getaa gens tis paler (40) 
E(¢) being the transform of the right-hand side of (39). 
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In the application of equation (40) it is necessary that the functions 
occurring should have a common region of regularity. We show that 
this is true for |n|<e where (=£+im and Imk=e. 

From the asymptotic form of H(«) we know that H(¢) is regular for 
n> but has simple poles at =+k. However, by inserting the factor 
(C2—k?), we may take (C2—k?)H(¢) to be regular for 7» 2—e. 

Moreover, by inserting the asymptotic form of I’ in equation (37) we 
obtain 


NGS eee 
Kq)~(z) =e On 


so that K(Z) is regular for | |<e. 
The explicit relation for K(¢) can be obtained from that of I(¢, y—y’), 
which is shown in Appendix I. to be 


MG y-y)= Pee can eee AS) 


where the argument of (42—2)! is taken to be 7/2 as | ¢|— oo on the real 
AXIS. 
It follows, from equation (37), that 


K(f) eb Gee CALL he ee 


where 


_ D—exp (2bi(—2)4}1 


Me 2bi(2—L) 


(44) 


and L(@) is regular for |y|<e. Moreover L(é) is non-zero inside the strip ; 
; . (nx? b : ! : 
its zeros occur at C= +7 ( = —k*}) and these points lie outside the strip. 


Thus the left-hand side of (40) is composed of two factors, L(¢) and 
b(k?—¢)H(¢), for which |7|<e is a common region of regularity. 

Finally we note, from (34) and (37) and the definition of E(x) contained 
in (39), that E(x) represents the y-derivative of the potential at y=b. 
Its asymptotic behaviour can therefore be derived from that of ¢,, outside 
the channel. With the help of (4) we obtain, for large 2, 

exp [tka] 
E(z)~—pofo) SPE on a le Soca aie ned 

_ so that E(¢) is regular for N<e. 

The next step in the solution of (39) is to identify L(Z) as the quotient 
of two functions, L,(¢) and L_(£), which are regular and non-zero for 
n> —€, n<e, respectively. Equation (40) will then become 


O(e— CH (C)L (C= ECL (2), - . > 3. (46) 
the left-hand side being regular for 7>—e and the right-hand side for 
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7<e. By analytic continuation they define, together, a function regular 
in the whole finite domain. It will be shown in Appendix IL. that the 
integral function so defined is a constant so that a knowledge of L.(Z) 
implies a knowledge of H(¢) apart from a constant factor. Reference 
to equations (26) and (29) shows that the expressions for the reflection 
coefficient R, and the function G(@) describing the transmitted energy, 
are independent of such a factor multiplying H(¢) so that its explicit 
evaluation is unnecessary. 

The functions L,(¢) and L_(¢) are obtained by applying Cauchy’s 
formula to log L(Z). Since L(Z) is regular and non-zero in the region 
In|<e, we can write, for |1|<e, 


log L(t) d 
log L(t) 5 SO eee (47) 


the integration being taken round an elongated rectangle within the strip 
and on the branch for which log 1=0 (fig. 3). 


Fig.. 3. 


From (44) we find that log L(t)~log|t| as |¢|— 0 within the strip. 
Consequently the contributions from the vertical sides of the rectangle 
tend to zero as the contour extends to infinity and we have the relation 


. 1 ¢ log Lid 1 log L(é) 
log L)= 5 | Ee eee 8) 
+ — 
=log L.,(¢)—log L_(2). Sree mete eee. (ST) 


The contours C, and C_ are shown in fig. 3, and, to ensure the individual 
convergence of L, and L_ we take 


fen log L(t) dt 


Aa, Pe Sel 
re ed a (28) 


| log L(t) dt = lim 


Cy t—{ T>o 


where 0<6<e. aa 

The second integral is defined in a similar manner. 

The two functions L,(¢) and L_(¢) so defined are regular and non-zero 
in the regions 7>—6 and 9 <6 respectively. In the limit when the two 
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; 
contours approach the axis, the two functions L, and L_ are represented! 
by the common expression 


Salle ven) 
sind al t—C 


eae ade i{1— exp ae (49) 
=en[ zlim | 8{ —anae—er SI 
and 


AP 4h = Lio 
For real k& (i. ¢. as «> 0), the contours C,, C_ should pass the branch 
point =k on the lower side of the real axis and t= —k on the upper side. 
The two contours then lie wholly in the respective regions of regularity 
of the two functions L,, L_ and the phase of (k?—#)* changes from 
0 to 7/2 for |t|<k, |t|>4 respectively. 


(50) 


Fig. 4. 


The above expression for L., could now be used to deduce the function 
H(¢) apart from a constant factor, and hence the asymptotic behaviour: 
of the reflected and transmitted waves explicitly. It is possible, however, 
to express L,(¢) in a somewhat more significant form which is most 
easily obtained by a deformation of the contour C,. The discussion of 
convergence is also simplified if we proceed from the derivative of log L_,(2) 
which is 


d SOLE i{1—exp [2bi(k?—#)*}}) dt ts 
EHO = safe, 8 — ae eage + 

Remembering that ¢ lies in the upper half plane, we deform the contour 
C, into the contour C, C,C, shown in fig. 4. The section C, is the 
upper and lower sides of the branch line extending from —0o to —k; 
C, consists of two ares of a circle, radius T, in the lower half-plane, while 
C, excludes the singularities of the integrand in (51 ) at the zeros of 
{l—exp [2b1(k*—?)*]} on the imaginary axis. 
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Along the upper and lower sides of the branch cut, the argument of 
(k°—t?)* is +7/2 and —z/2 respectively, so that the contribution to the: 
right-hand side of equation (51) from C, is 


Els ia a [wes eal dit 


Ini 2b(k2—B)) ¢—b: 
sa 19 s|= 1—exp [—2b(?— “aH | dt. BE f* (P?—k*)8 dé 
~ Bai (I—exp [—20(F—F)) Sigel (aaa 
gos b = (7? — j2 ee) cosh? T/k du ts 
= =| Spe + cosh"? T/k— ae EeESR EE: poarat| - . (52) 
so that 


T b 2¢ De ee \a “ 
anae el + = cosh Z| =| i+ gern (a) |: . (53), 


Consider now the integral round C, in the quadrant for which. 
—m<argt<—7n/2. In this range —7/2<arg (k?—#)}<0, and the 
asymptotic value of the integrand in (51) is —(log 2bt)/t? so that, in the: 
limit T + oo, the integral along C, is zero in this quadrant. 

In the range —7/2<t<0, however, we have 0 <arg (k?—??)?<7/2. 
On this section of C, the asymptotic value of the integrand is 
(log exp [—2bt])/t?—=—26/t and the contribution is therefore 


0 
-- (¢- =| Re meme a = (Bay 


Oni TV —n/2 2 


Finally, there is the contribution from C,; due to the phase difference of 
the integrand on the two sides of the imaginary axis. Since the zeros: 
of [l1—exp {2bi(k?—#?)*\] occur at t= —i(n*x?/b?—k*)?, it follows that the- 
integrand in (51) has the same phase, along that part of the negative real. 
axis for which ear. as 


where N is the largest integer not exceeding Tb/z. It foliows that the- 
integral along the two sides of the negative real axis is 


4 ace [2b2( oo Le dt 


Oni) ¢ a (iO? 


=— ee it 1 {e4i(5 


= n=1 


—iT dt N iL 1 ‘ 
3) ps ebay ah 2 fia lak (neqee.. VEE 
n=1 -i( a -i) #(¢—¢)? wan bal +i("p —i#) 
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If the contributions from the various sections of contour, given by (53), 
a (55), are now added, we obtain 


as LiKe p =f 26 ans Ue J 3 
7 i | enna) CEI ra 5 
b ed 
oh y + = cosh zh. 


N 
+2 lim [ 2 dian Se eg en 
T+ RY oe c+i( 2a? —i#) 


Now, since, for large z and N, 
cosh-1z~log 2z 
N 


1 
a = ~log N+y, 


n=1 


where y is Euler’s constant, it follows that (56) can be written 


_, (k—-“\ b 
<rlog Ly, =2| 4+ ga 7 tan (a) | airs 


N 1 tT b 
+ lim| 2 (ova amon 
Tivied | wel joa (F -is)) C+i1T am 
ee ea Tie bone 30 fhe ae 
= 1 b b : ; 
eee |\—a=— ae +24 log, = rol ek) 
n=1 Ean 3 —i*) 7 7 a7T 


since N~Tob/z. 
A simple integration with respect to ¢ then gives 


L,W)=Kexp| — $+ “ogre — Hiv, MP 2 w—eytean (FS) 


the k2b2\47] — ibe 
+ B tose +4( at) lapel | 
payghe k2b2\4 Lb itb Lb ae b 
-« f (E8) 2] (2) on] 2 (otc) 
+ ee eyitan (73). fear LTH =! 1345) 


Since 


equation (58) gives 


L (LQ) =L)=(«'y IT | 1 | exp (002) 


(«")?6 


= ee —exp [2bi(h?—@?)#]), 5 2 . (59) 
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nd comparison of equation (59) with (44) shows that the integration 
constant «’=1. 
Moreover 


kb2\) ig igo 
ees HR GCL on cee NT 
1 = b\t 10d ith 
es aft : e{(1 — aa) - = 2 =| 


cade z lo g{i— Ln a3 eli > ey 
a= n=1 re) nT 


B 
so that i. 
= ee | ex ae (1-1 Hog ii) 
ce ae (2—@)3 tan-2 (=) ad) 5 ea 3 zl , 
so ( 2 -i#) | 


and 


sin (bk sin 6) ]! bkcosé | ibk cos6 /, Qa 
..(& cos ¢)= | SES | “| — ee i (1-y-+108 77) 


(bk 6 si 20 k 6 bk cos 6 
ibkOsing plete iaLt ee By NEOSe | a 
n=1 cs LSB < gun ae 
a) { 
ibk Qa eee ewe OK OK 
(i)=exp| —F += (1 —y-+log zz) —i © {sin fu hee 


It is now a simple matter to obtain expressions describing the reflected 
nd transmitted waves. The reflection coefficient is, from equation (26), 


Res ResHit)} Res HC). ae Ow ef. (64) 


=-k c=k 
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and since, (see Appendix) 
(k= 22) (OU (0) Dyer ok. ee 


where D is some constant, it follows that 


Res =+ v 
ves = Ta TAC 
c=4+k 2kL4(+4) 


Equation (64) can now be written 


R=—|R|exp [2ék]=—L,(b)/L,(—k)=—{L,(k)}2, . (67) 


where / is the end correction. 
Using the expression for L(k) given in equation (63), we find that 


(BR l=oxp [- bE) one et eee 
1 2a 1 eae OK 
Baviad bt ck ams} ee 1 ait Bea . 2 (6SF 
=-(1 ree in) hee {sin oes Sd 
Finally we derive the expression for G(#). Equation (65) shows that 
D 
————————. . . . . . (70 
rE (k sin 6)?L, (& cos @) a: 


Substituting (66) and (70) in (29) we obtain 


2 sin? (bk sin @) | La(k) P 


C= bksin?@ = =(1—| R|)?| Ly(kcos@) > © at 


__ exp [bk cos 8] sin (6k sin 6) (72) 
at sin 6 sinh bk : Se 
by virtue of relations (62), (63) and (68). 


A simple check on this last relation is obtained by noting that from its. 
definition the function G(#) must satisfy 


i Gi) ofan hes OOF ieee 
To verify this we consider the integral 

5 fexp [bkz] f 

: gees dz, . . - . . . . (74) 


taken round the unit circle in the z-plane, with indentations at the poles. 
Z2=+1. 

The integrand is free from singularities inside the contour and the 
integral is therefore zero. If we substitute z=exp [76], then this is. 
equivalent to 


Or 
_— 


P | exp [bk exp (10)} sa —Ti | Res + ne} = (a eee 
J —n7 - z=-1 
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Or 
iy sin (bk si 
| exp [bk cos 6] Se Oren) 8a (exp [bk]+exp [—bdk]) (76) 
so that 
E CAG an : sin (bk sin @) 7, | 5 
J _4@) d= | bee exp [bk COs 6] sin @ sinh bo ane (77) 


DISCUSSION OF RESULTS. 


One of the most striking features is the extremely simple relation for 
‘the amplitude of the reflected wave (equation (68)). The corresponding 


0 02 O4 agi 06 08 -0 


: d bk 26 
The variation in end-correction as a function of eAU aoa 


result for the circular pipe is more complex, but, for large values of the 
wavelength, the amplitude is approximately proportional to exp [—k’c?/2] 
where c is the radius of the pipe (Levine and Schwinger 1948). The 
decrease amplitude with the wavelength, when the latter is large, is 
therefore smaller than in the two-dimensional case. 

The end correction is illustrated graphically in fig. 5. Hquation (69) 
shows that it increases indefinitely as the wavelength becomes infinitely 
great, as opposed to the circular pipe for which it tends to a finite limit. 
The same conclusion was reached by Rayleigh (1904) using approximate 
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methods. The physical significance of the end-correction is consequently 
somewhat remote and it is more constructive to consider the actual 
correction to the phase in the returning wave which is given by 


2bk Qa 2 bk bk) 
Sblan ne (Tes i bo sin-t-— =), ee 
2k e (1 y+log in) 2 E {sin ai _ (78) 


and tends to zero for large values of the wavelength. 


Fig. 6. 


0 0:2 04 pp 06 08 1-0 
The radiated energy as a function of 6/7 for various values of bk/z. 
In fig. 6, the variation in G(@) is shown for various values of bk/a( = 2b/A) 
When the wavelength is comparable with the width of the channel, most 
of the energy is transmitted in a forward direction. The function 


oo 
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decreases monotonically with 6 except for values of bk/m near unity when 
there is a minimum in the region of 6=7/2. If bk—=r, for example, 
G(z/2)=0 and no energy is propagated in this direction. Beyond this 
value G(@) increases slowly to 0-0118. 


APPENDIX I. 
To derive the Fourier transform of 


i : . 
Ea; y)= 5 HP ha t+y*)"] Se eaenuiett. be GAS) 
we begin from the known relation (Watson 1944) 
HWY@)= =| exp [iz cosh ¢] dt: if O<argz<n, . . (A. 2) 
7 —"00 


so that 


Tie: y=1[ exp [tk(a?-++-y?)? cosh t] dé if Imk=e>0 


I 


o—10 
4 ( exp [kx cosh u—k| y|sinhu]du, . (A. 3) 
J —w io 


where tan =| y|/z. 


Figy 1: 


Substituting ¢=k coshw this becomes 


1 pexp[eal+il y |(e—&) dl 
| eae ei gee. (AG 4) 


where the contour of integration begins and ends at infinity and passes 
between the two branch points C¢=-+k (see fig. 7). By a deformation 
of the contour, this can be written, if |7|<e 


(a, y)=— oh Bethe (k2— 4) 


(A. 5) 
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It follows that the transform of l(a, y) is I’(¢, y) where 


apse Gare Hs SAR 


and is regular for | Im |<e. 


APPENDIX II. 


During the course of the analysis, use was made of the fact that the 
integral function defined by the two sides of the equation 


b(k*—@)*H(C)L ,(C) = E(C)L_(¢) 


was, in fact, a constant. 

To prove this we examine the asymptotic values of the two sides, in 
their respective regions of regularity, as | ¢|— oo. 

For the asymptotic behaviour of L,(¢) we refer to equation (51), 


ee ih tt T 1{1—exp [2bi(k?—#*)*}}) dt 
pool sain | 8 area |=] 


Ree. Weegee i{1—exp [2bi(k?—#)*]}}) dt 
=exp Ea: ae ee \< 


£7?) flexp [—2b(@—)4}|_ dt 
Now 
ft |, ffl cexp oi(e—e yp 
aa mice See IS =—0 (3) . (A. 8) 
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log (== eee Pp did 
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3 (A. 9) 
my i 0 1 
Substituting in (A. 7) we obtain 
i 
L,(l)~— (A. 11) 


and 


Se (ID) eee, (A. 12) 


oe 
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For the asymptotic behaviour of 


HC) =f H(x) exp[—ilx]dz Imf>0, . . (A.13) 
cand 


} E(g)= i” K(x) exp [—ila] du Imf<0, . .-. (A. 14) 


‘we note that the significant contributions to the integrals, when | ¢| is 
large, will arise from the values of H(x), E(x) for small negative and 
‘positive « respectively. Now the difference in potential across the 
plane y=b must be a continuous function of « which vanishes for «>0. 
It follows that, as «> —0 


eerie a? Os ee os wee tN 1D) 

or 
Time Giece ter eet te. GAS 16) 
Again the y-derivative of the potential across the plane y=6b must have 


_.an integrable singularity at c=0, and since this is just E(x”), «>0, we have, 
as «> +0 


L(G ear eeo eall sperm emi Se (Ad) 
“Or 
E(Q)~(i0)P-1. teres... (A..18) 


Using (A. 11), (A. 12), (A. 13) and (A. 14) we see that the behaviour of 
the integral function is, in the upper half plane 


b(8—PYH(OL, (Qi), 
while in the lower half plane 
B(Q)L_(Q)~ (io). 


Because of the limits on « and £, the integral function cannot tend to 
infinity faster than ¢* and must therefore be a constant. 
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ABSTRACT. 


The magnetic design of betatron-starting flux bars for electron 
synchrotrons is considered in relation to the principles governing efficient 
initial electron injection and efficient betatron-synchrotron transition. 
A brief discussion of betatron-synchrotron transition theory is included, 
in order to show how flux bar performance may be specified therefrom.. 
A design procedure is derived for determining the essential flux bar 
dimensions, approximate corrections being made for skin effect in the 
laminated material of the bars. 

To illustrate the principles discussed in the main paper, the results 
of some detailed performance calculations on the 375 MeV. Glasgow. 
University synchrotron are given in an Appendix. 


$1. INTRODUCTION. 


In designing the magnet of a betatron-started electron synchrotron,. 
economy in weight and space dictates the use of betatron flux bars of the 
minimum size consistent with satisfactory performance. LHssentially 
this implies that the betatron acceleration phase of the machine is 
restricted to imparting just sufficient electron energy to ensure efficient 
transition to synchrotron acceleration. The betatron bars may then 
saturate, while synchrotron acceleration maintains the energy increase. 

Two fundamental criteria govern satisfactory betatron bar performance. 
The first is that, to obtain efficient electron injection into the betatron 
phase, the bars must produce an equilibrium orbit located centrally in 
the vacuum donut. The second is that, during the betatron phase, the 
electron orbits must be maintained sufficiently central in the donut to 
permit efficient transition to synchrotron acceleration when the R.F. 
accelerating voltage is applied. As Goward (1949) has shown, this. 
second criterion involves marked differences between machines of various. 
energies. 

For the sake of definiteness, attention in this paper is focused primarily 
on synchrotrons producing from about 100 to 400 MeV., and having a 
general magnet arrangement of the type illustrated in fig. 1. With the 
exception of §§3 and 4, however, the treatment throughout applies, 


* Communicated by Denis Taylor. 
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without essential modification, to other energy ranges and magnet designs. 
In §2 it is shown how the first performance criterion mentioned above 
determines the dimensions of the air gap (or gaps) controlling the betatron 
flux. The second criterion, namely efficient betatron-synchrotron 
transition, effectively determines the minimum permissible betatron bar 
cross section. §§3 and 4 show how the simple transition theory given 


General arrangement of a synchrotron magnet. 
. Air gaps controlling betatron. flux. 
. Betatron flux bar. 
Main excitation. coils. 
. Vacuum donut. 
Magnet poles. 
. Magnet return path. 


KHOOW 


by Goward (1949) may be put into a convenient form for specifying 
betatron bar performance, treating the problem in full only for synehro- 
trons of the energy range mentioned above (100 to 400 MeV.).  §§5 and 6 


D2 
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then give a method for determining an appropriate betatron bar cross 
section, employing the D.C. hysteresis loop of the bar laminations, and 
§7 shows how a moderate amount of “skin effect” in the laminations 
may be allowed for. No account has been attempted of the effects of 
cross fluxes in any part of the magnetic circuits, although in practice 
such fluxes are difficult to avoid entirely. 

In Appendix I electrical methods of modifying the betatron-starting 
flux in a synchrotron are briefly reviewed. In Appendix II, some 
detailed performance characteristics of the 375 MeV. Glasgow University 
synchrotron are presented, illustrating various topics discussed in the 
main paper. A particular attempt is made to estimate overall electron 
trapping efficiencies more accurately than is possible by Goward’s much 
simplified theory. 

§2. Betatron Arr Gap DESIGN. 
(i) Fundamental Equilibrium Orbit Equations. 

To establish an equilibrium orbit of radius 7» in a betatron, the well 

known Widerée flux condition must be satisfied : 


Odg/ Ot 2rr20H [OF ey he: anos 1.) es 


H, is the vertical field at radius ry in the annular magnet gap, and 
dy is the total flux linking the orbit. Internal leakage fluxes of the 
synchrotron magnet, together with the flux in the annular gap at radii 
less than 79, contribute an appreciable fraction of ¢9, and may be expressed 
as (1—k,)2ar2H», where, for a given magnet, ky depends only on 7). The 
remainder of ¢, is the flux ¢, carried by the betatron bars. Thus 4, and 
ry are related by the equation 


Obs Ot—2arak GH y/Ok | Le. & Sas ae 
Assuming that 7) denotes some fixed radius, a slight adjustment of the 
bars will change 0¢,/0t in relation to 0H)/dt, and ¢, will satisfy an equation 


Ob, /Ot=Qargky(1-tn)OH pO. . . . . . . (8) 


The equilibrium orbit will now be of radius 7)+6r9, where, provided 
y and (679/79) are both small, 

Ori to kag (ha) See oe eee (4) 
The field in the annular gap is assumed to follow the usual r~” radial law. 

To illustrate the orders of magnitude involved in equation (4), values 
appropriate to the 375 MeV. Glasgow synchrotron are n=0-7, 77> =125 em. 
ky=0-82. It follows that a misadjustment in the rate of change of ae 
in this machine of only 0-3 per cent involves an error in Squihhonan 
orbit location of 1 em. 

(ii) Magnetic Circuit Equations. 

To express the Widerée flux condition in practical terms we must con- 
sider the physical layout of the synchrotron magnet, e.g. as shown in 
fig. 1. For simplicity of analysis any actual arrangement of betatron 
bars may be represented nearly enough by a single laminated bar of 
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length L, useful cross section A,, and carrying a uniform flux density B,, 
placed in series with an air gap of length g, area A, in which the flux 
density is B,. Thus for continuity we have 


6 )=A,B,=A,B,. Say TRIES ores gar nek CD) 


Let G be the height of the main annular air gap of the synchrotron, 
at radius 79, and let H,, denote the magnetic intensity in the betatron bar. 
Assuming equal return path reluctance for annular and betatron bar 
fluxes we have , 


PereeuH, =CH ee Dee et te no (8) 


where y’ =gA,/A, and the term D represents any D.C. bias which may be 
applied to the betatron bars. 

When skin effect in the bar laminations is negligible, B, and H, are 
simply related by a normal D.C. hysteresis loop taken well into saturation. 
When skin effect is appreciable, B, becomes the average flux density in 
a lamination and H, denotes the intensity at a lamination surface, and 
some modification of the hysteresis loop is required to give the correct 
relation. This problem is not considered further here, since the results. 
of §2 (ili) are not affected, but a method of constructing a suitably 
modified loop is given later in §7. 

Differentiating equation (6) partially with respect to time, and using 
equation (5) gives 

hepa teins = Ke yd xy al Gk A Sey mm 
where R is defined by : 
BG RUS ee ee cs fh ee (8) 


and s is the “differential permeability’ 0B,/0H,. It follows from 
equation (2) that an equilibrium orbit of radius 7 will be established when 


R satisfies the equation 
G/R=2zrzk. SE Hes Natta (9) 


(iii) Betatron Air Gap Dimensions. 
To determine the betatron air gap dimensions it is permissible to assume 
that, in equation (8), (L/sA,)<g/A, at low flux densities. Equation (9) 


then approximates to 
g/Ag=G/ (27g). PES oay «a une 6 (LO) 


For economy g and A, are made as small as possible, the limiting factor 
being the tolerance on the dimension g. Equations (4) and (7) show 
that a small error, 5g, in g involves an error in equilibrium orbit location 
given. by 
Of 9/9 = —Kdg/(1—n)g- MeeeM ay en wuhatia (ll) 
Observations on a 30 MeV. synchrotron (Fry et al., to be published) 
suggest that, if ry denotes the mean radius of the vacuum donut, within 
90 per cent of the maximum electron injection into the betatron phase 
will be obtainable provided 67, is less than 0-1 of the radial donut width. 
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Applying this information, for example, to the 375 MeV. synchrotron, 
which has g=1-9 cm. and a donut width of about 12-5 cm., 679 must be 
less than 1:25cem. From equation (11), g must be correct to 0-007 cm. 

Thus equations (10) and (11), plus the requirement for accurate location 
of the equilibrium orbit at injection, decide the values of g and A,, and 
the value of L will generally follow geometrically from these. So far the 
discussion has involved only a slight restriction on A,, namely that the 
approximate equation (10) should be justified. To specify A, more 
definitely requires some further specification of betatron bar performance, 
i.e. the rate at which the electron orbits may contract, as a result of 
bar saturation, while still giving efficient betatron-synchrotron transition. 


§3. BETATRON-SYNCHROTRON TRANSITION. 


(i) Transition Theory. 

The problem of electron trapping at betatron-synchrotron transition 
has been discussed in several papers on synchrotron theory (e.g. Bohm 
and Foldy 1946), and particular aspects, such as the effect of finite time 
of rise of the synchrotron R.F. voltage pulse, have been studied (DePackh 
and Birnbaum 1948; Twiss and Frank 1949). A simple treatment 
which is sufficiently comprehensive for the purposes of this paper is that 
given by Goward (1949), who assumes sudden cessation of betatron 
acceleration and simultaneous application of a steady R.F. voltage. 
Goward shows that quite different factors determine efficiency of 
transition in a synchrotron giving say 30 MeV. and one giving a few 
hundred MeV. ‘To avoid excessive length, the discussion in this section 
(§3 (i)) is confined to synchrotrons giving from about 100 to 400 MeV., 
other energy ranges receiving only brief consideration in §3 (ii). 

Following Goward’s treatment, we assume that electrons are accelerated 
by betatron action up to a total energy E,. An R.F. voltage of peak 
value V is then switched on, trapping a fraction f, of the electrons into 
stable synchrotron orbits. In a high energy machine, Goward shows 
that V is always large enough to make f, approximately unity; 7. e. 
nearly all the electrons are trapped. The ‘trapped electrons execute 
radial oscillations about the equilibrium synchronous orbit (Bohm and 
Foldy 1946), ranging in amplitude up to a stable limit p,, where p, is a 
function of V. If V is increased beyond a certain value, V 1> Ps exceeds 
the clearance (p,,) between the equilibrium synchronous orbit and the 
donut wall, and some trapped electrons are consequently lost. V, and 
P» are related by the equation 


Vi =[1—n)PirK Bonito), (12) 
where K is given by 
K=1-+][(1=n)pe Pe ee Sees) 


8, is the electron velocity, at total energy E,, relative to that of light. 
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For V>V, a fraction (1—f,) of the electrons collide with the donut walls, 
while f, are accelerated further ; f, is given by the equation 


p= (Birr RIF EAN 5/ Mil Page Ps tsiotaiee me le re 4) 
shown graphically in fig. 2. 
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Electron trapping efficiency factors. 
fo. Result of scatter to donut walls at transition. 
fz. Result of change of stable R.F. phase during 
synchrotron acceleration. 


As synchrotron acceleration proceeds the betatron flux bars saturate 
* and the accelerating ‘“‘ volts per turn ”’,.v, which must be supplied by the 
R.F. resonator, rise to a maximum vp). Assuming that electron radiation 
losses are not of primary importance in synchrotrons of the energy . 
range considered, vy is given by the equation 


p= 4refrHky|(3X10%),  . 9.9. 2... (18) 


where f is the synchrotron magnet excitation frequency in ¢/s. and E the 
final (total) electron energy. Thus V>v is a necessary condition for 
synchrotron acceleration and Goward introduces a further efficiency 
factor fs, depending on the ratio V/vy, which takes into account electrons 
Jost during synchrotron acceleration. The factor f; is also shown in fig. 2 
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and is calculated for high energy machines by neglecting phase damping: 
of the synchrotron oscillations (Bohm and Foldy 1946), This damping is,. 
of course, important, but is difficult to include in a sufficiently general 
manner. ‘To a sufficient approximation, the fraction of original betatron 
electrons accelerated to peak energy is given by f. or jf, whichever is the 
least, and hence the most efficient value of V to use is one which makes. 
fo=fs, e.g. for 65 per cent efficiency, putting f,—f,—0-65, fig. 2 gives. 
(V/v,)=3°3 and V/V,=1-35. Thus we have 


ViG=2-45 0 oe O.., SO 


By means of an example, it is shown in the next section that an equation 
such as (16) forms a convenient basis for specifying betatron flux bar 
performance. It must be emphasized, however, that the equation is. 
based on essentially pessimistic calculations of electron trapping efficiency... 
Even a simple rectangular R.F. pulse has in practice an appreciable 
time of rise, and phase oscillation damping will reduce loss of electrons 
to the walls. The same damping appreciably reduces loss during’ 
synchrotron acceleration, and thus both f, and f,; tend to. become nearer 
unity. Moreover it is possible to reduce both kinds of loss by suitable 
shaping of the R.F. pulse (DePackh and Birnbaum 1948; Twiss and 
Frank 1949). <A relatively small R.F. voltage would suffice for initial 
trapping of the electrons (f,—1), after which, to allow for the increasing 
volts per turn v, the applied R.F. voltage would have to be increased.. 
This latter process in itself gives increased oscillation damping. 

(ii) Flua Bar Performance. 

Taking the 375 MeV. synchrotron as example, the parameters in 
equation (15) are f=50 c/s., r>=125 cm., E=375 MeV., k,=0:82. Hence 
¥=2-5kV., and from equation (16) the condition for 65 per cent efficient 
overall electron trapping is V,=6-1kV. In equation (12), however, 
n—0-7, and assuming 8,1 gives 


Vis (S/L0P nee Guts lee, class ee ee 
Since V,=6-1 kV. we have 
Ep? = 203108. 1-05. eee ee ee 


Equation (18) defines pairs of values of E, and p,, which permit 
65 per cent efficient transition, on the basis of the above simple theory.. 
Since the effective donut width is 12-5 cm., if we write (6-25—p,,)= Ar, 
the distance dr represents the permissible radial displacement of the- 
betatron beam centre from 79, at transition energy E, (assuming that 
transition occurs when betatron beam and stable synchronous orbit 
coincide). A table of values of E, and 4dr may now be drawn up, as in: 
Table I. 

TABLE I. 
E,(MeV.) 5:2 6-2 7:4 9-0) 11:3 
Ar(cm.) 0 0-5 1-0: 15: 2-0 
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Thus in addition to establishing an equilibrium orbit of radius 7,, when. 
initial electron injection into the betatron occurs, the betatron bars must 
satisfy the additional requirement of maintaining the electron orbits. 
constant within the limits of Table I. The method of extending the bar 


_ design to meet such a requirement is described in the following sections.. 


It is convenient at this point to consider briefly how the discussion in 
this section (§3) must be modified for synchrotrons of other energies.. 
Above 400 MeV., electron radiation loss becomes important, and may 
determine the maximum R.F. “ volts per turn” v) (see Appendix II).. 
Moreover R.F. pulse shaping becomes essential to avoid low f, values. 
(DePackh and Birnbaum 1948) and the values of E, and Ar will be governed 
by a more complex transition theory than that of Goward (1949). 

Considering energies below 100 MeV., the values of v) and of V become 
smaller, and loss to the donut walls at transition (f,) becomes less likely. 
The value of Ar is not so critical, although Ar must always be less than 
half the donut width (7. e. the stable synchronous orbit and betatron 
beam must always coincide within the donut limits). The factor /,,. 
however, which depends on E, and V, now becomes more important than 
J> or f3, and the design problem consists essentially of arranging a high 
enough value of E, to make f,;—~1, taking into account the value of V 
available. A full treatment is given in Goward’s paper. 


§4. Berarron Firux Bar Cross SEctTION.. 
(i) Instantaneous Circle Orbit Equations. 

Neglecting free oscillations, individual electrons in a betatron beam. 
move on so-called *‘ instantaneous circle ” orbits (Kerst and Serber 1941) 
which, at any time, have a range of radii corresponding to the energy 
spread of the electrons. It is convenient, for the following discussion, 
to consider the electron motion in terms of these instantaneous circle 
orbits, and to use the instantaneous circle orbit equation (Amaldi and 
Ferretti 1946) : 

6,= 2nr"H,-+C. 5 Ree Ge, aie AMD 


H, is the vertical field at radius r and ¢, is the total flux linking the 
circle of radius 7. If C is a constant, the equation defines the variation 
of r for any electron throughout its betatron motion. 

For design purposes it is sufficient to consider an electron injected from 
the electron gun so that, initially, r=r). (This electron roughly defines 
the mean orbit radius of the whole electron beam.) At injection, equation 
(19) for this electron takes the form 
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After some betatron acceleration, 7 will in general have changed to some 
value 7)+4r, where we assume (47/r9) is small. Corresponding to this,. 
the flux through the circle r will then be given by 


ibe Oar He ee ete eee ee (21) 
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where (A¢,/f>) is also small. 4¢) and Ar are related by the equation 
Ado|(27r2Hy)=(1—n)(Ar/ro). =. - + + + (22) 


(ii) Determination of Betatron Bar Cross Section. 

From §2 the betatron air gap dimension and bar length L (as defined 
in §2) are already determined. For the present purpose, Hy is nearly 
enough zero when injection into the betatron occurs. Thus when Ho=9, 
let the variables in equation (6) have the values B,=B,, H,=H,. 
‘Then for any value of Hy we have 


B,—B,=(G/g')Ho—(L/g’)(H,—Hy). « - + (23) 


Equation (23) gives the total change in B, corresponding to an increase 
in field at 7) from zero to Hp. 

A theoretically ideal betatron bar would have constant differential 
permeability s, and would maintain the electron considered in §4(i) at a 
constant orbit radius r=ry. Let the parameter R in equation (7) have 
the value R, when B,=B, (and s=s,). Then, for an ideal bar, R would 
be constant at R,, and equation (7) gives 


B,—B,=GH,/(A,R;). . . . . . . . (24) 


Comparison of (23) and (24) gives the amount, 4B,, by which the actual 
bar density differs from the ideal value, as a result of variation of s and R. 
Assuming all non-linear effects to be confined to the betatron bar flux, 
we may identify A,4B, with the 4d, of equation (22), and obtain 


(Ry/L) 4¢y=(B,—B,)(1/s,)—(H,—H,). - + + (25) 


If we neglect skin effect for the moment, the right hand side of equation (25) 
has a simple graphical meaning in terms of the D.C hysteresis loop as 
shown in fig. 3, and may conveniently be denoted by (—4H,). Using 
equations (9) and (22), we have 


— 4H ,=[(1—n) G/L Hy dr 9 eee) 
Further, using the well known relation 
EH, B= 300rgHy, ies te ty cape ee) 
equation (26) may alternatively be written 
— 4H ,=[(1—)G/(300 Lr2k,)JE,B,47. . . . . (28) 


In conjunction with a pair of values of E, and Ar, specified for example 
as in §3, equation (28) effectively determines a betatron bar cross section. 
The derived value of 4H, specifies pairs of values of B, and B, from the 
hysteresis loop of the bar material, as in fig. 3. From equations (9), 
(24) and (27), A, is then given in terms of such a pair of values by the 
equation i 


A, =2nr Hh ,8,k9/800(B,—B,). ... . . ~ (29) 
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Magnetic Flux Density (B) 
9O 


Magnetic Intensity (H) 
(a) Little (or no) Negative Flux Bias. 


Magnetic Flux Density (B) 


A, ; se 
Magnetic Intensity (H) 


(6) Large Negative Flux Bias. 


Hysteresis loop construction for determining changes in electron orbit radius. 


In the following section, the application of the above procedure is 
illustrated by a specific example. 
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§5. EXaMpLe DesiGN : 375 MEV. SYNCHROTRON. 


(i) No Applied D.C. Flux Bias. 

Taking again the example of the 375 MeV. Glasgow synchrotron, the- 
known parameters in equation (28) are n=0-7, G=10 cm., L=100 cm., 
ry =125 cm. and ky=0-82. If at least 65 per cent efficient transition 1s 
taken as a design criterion, Table I. gives appropriate values of E, and 4r.. 
For example, the values E,=9-0 MeV., 4dr=1-5 cm., inserted in equation 
(28), give —4H,—0-105. Drawing a tangent at the point B,;=0 to the 
hysteresis loop of the transformer steel laminations to be used yields the 
value B,—4600 gauss, in the manner indicated in fig. 3 (a). From. 
equation (29), A,=4200 sq. cm. - 

Other values of E, and Ar from Table I. give cross sections greater 
than 4200 sq. cm. Hence we conclude that if no D.C. bias is to be used, 
and skin effect in the laminations is negligible, a cross section of 4200 sq. em. 
is the minimum which will ensure at least 65 per cent transition efficiency. 
(ii) Application of D.C. Bias. 

An auxiliary D.C. winding on the betatron bars may be used to give 
B, some value other than zero, and to effect a saving in bar cross section. 
The graphical construction may now take the form of either 
fig. 3 (a) or 3 (6), according to the amount of bias used. 2000 gauss. 
negative bias, for example, gives a fig. 3 (a) type of construction, and the 
same design procedure as in §5 (i) yields a cross section of 3200 sq. cm. 
A greater bias gives the type of construction shown in fig. 3(b). Here, 
as is readily seen from equation (26) or (28), the electron orbit, 
initially of radius ro, will first expand to some extent, and then contract 
through ry as tangent and hysteresis loop intersect. The design procedure 
is unaltered, and taking B, as 4000 gauss negative gives a cross section 
of 2000 sq. cm., less than half that determined for zero bias. 

Equation (26) may be used to check that the initial orbit expansion 
under fig. 3 (b) conditions is not excessive ; e.g. in the above 4000 gauss. 
bias example, instantaneous circle orbit expansions from 125 em. at 
injection to a maximum of 126 cm. at about 3-6 MeV. electron (total) 
energy will occur. 

By the above methods a betatron bar cross section can be determined. 
to give any specific degree of electron orbit contraction. One single 
value of cross section is obtained if no flux bias is to be used, but a variety 
of combinations of cross section and bias will give essentially the same 
results. So far the question of skin effect in the bar laminations has. 
been neglected and we now proceed to determine whether practical 
lamination thicknesses involve any serious modifications in betatron: 
bar performance. 


§6. Skin Erect, 
(i) Modification of Hysteresis Loop. 
The use of practical lamination thicknesses in betatron bar construction: 
will in general involve appreciable skin effect, and the variables B, and H,, 
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used in §§4 and 5, should strictly be defined respectively as the average 
flux density in a betatron bar lamination and the intensity at a lamination 
surface. Some modified hysteresis loop is therefore required to relate 
B, and H, if the graphical method of §4 is to yieid a correct result. 

Consider a lamination extending in thickness between the limits 
—§/2<«<€/2. The flux density, B, and intensity, H, at any point in 
the lamination (directed say along the y axis) are considered to be 
functions of x and time, t, only, and are related by the well known 
equation 


@?H/da2=(47/10%)OB/At, . . . . . . (30) 


where p is the electrical resistivity of the material. Solutions of the 
equation when B is a linear or sinusoidal function of ¢ are also well known, 
but to a first approximation the particular function of time is immaterial, 
and it is found that 0B/dt may be treated as a constant with respect to a. 
Equation (30) is then simply integrable, and it follows that the distribution 
of H is of the form (a+62?). If H, and H, are the values of H when 
“x=+€/2 and 0 respectively, then 


H,—H,=(a/10°.2p)20B,/at, . . . . . (31) 


where B, is the value of B averaged over x. 

Assuming that the synchrotron magnet excitation is derived from a 
resonant circuit of frequency jf, the magnetic field Hy will not vary 
sinusoidally, since saturable betatron bars introduce non-linearity into 
the circuit. If the ratios of magnet inductances with unsaturated and 
saturated bars is g, then during the initial betatron phase 


CHO aoe (82) 


where H,,, is the peak value of H,. Using equations (7), (9), (27) and (31) 


we obtain 
H,—H,=(2:07/10)rokyfE@/GAyp, . - - + (33) 
where E is the peak electron energy (eV.). 
Corresponding to the range of variation of H given by equation (33), 
the flux density B will vary from B, at the lamination surface to B, at 
its centre, where (H,, B,) and (H,, B,) are points on the normal D.C. 
hysteresis loop. If we approximate to this range of the hysteresis loop 
by a straight line, the distribution of B will also be of the form (a-+-62*), 
and by simple integration 
Bj Di (Pee octhemtnuemnsy, «os ~ (84) 


As indicated in fig. 4 this approximate result gives a graphical method for 
deriving from the D.C. hysteresis loop an (H,, B,) curve representing 
true operating conditions. 
(ii) Example: 375 MeV. Synchrotron. 

Again taking a specific example, the 375 MeV. synchrotron, the 
proposed bar material consists of 0-007 in. laminations of resistivity 
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50 microhms per cm. cube. With 7y9=125 cm., kj=0-82, f=50 c/s. and 
q=3, equation (33) gives 

H,—H,=830/A,; . . . . . . . (35): 


if A, is expressed in sq. cm. By neglecting skin effect, a cross section 
of 4200 sq. cm. was determined as suitable for use with no D.C. bias. 


Fig. 4. 
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Modification of hysteresis loop by skin effect. 
(a) D.C. hysteresis loop, 
(b) Modified loop By —B,=(Bs—B,)/8. 


For this cross section (H_,—H_,)=0-2, which is small enough for the approxi- 
mate methods of the last section to be applied. Deriving the modified 
loop, as in fig. 4, and applying the construction of fig. 3 (a), with B,=0, 
indicates orbit contractions slightly less than those predicted by the 
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unmodified loop. Neglecting skin effect resulted in an over estimation of 
the cross section required by only a few per cent. 

Similar eee may be made on the assumption of A, values. 
determined in § 5( The practical result is to predict that, for a given 
cross section, Eke less bias would be required than that specified in 
§5(ii). For example a 2000 sq. cm. cross section would require about 
3000 rather than 4000 gauss negative bias. The detailed variations in 
electron orbit radius are of course affected both by the change of bias and 
the modified loop shape, and rechecking for excessive initial expansion 
as in §5(ii), shows expansion up to 127cm. at 5 MeV. total electron 
energy. 

It is Puovk ble to conclude that in general, provided equation (33), 
indicates no excessive skin effect, use of a normal D.C. hysteresis loop is. 
sufficiently accurate for initial en purposes. <A final design should,, 
however, be checked by the appropriate modified loop. 


§7. CONCLUSION. 

The essential steps in the magnetic design of betatron bars for an electron 
synchrotron may now be summarized as follows. Betatron air gap. 
dimensions are determined from equations (10) and (11), and the flux 
bar length generally follows simply from geometrical considerations.. 
The methods of §3 define the rate of radial shrinkage of the electron 
beam orbit which may be permitted if efficient betatron-synchrotron 
transition is to be ensured. Equation (28) essentially expresses orbit 
shrinkage rate in terms of non-linearity of the betatron bar material, 
and a graphical construction on the hysteresis loop, as shown in fig. 3,. 
together with equation (29) gives the appropriate cross sectional area of 
bar. As a check on the validity of neglecting skin effects equation (33) 
should be evaluated, and a final design should be checked by proceeding 
to construct a modified hysteresis loop by the method of §6(i) and fig. 4 
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APPENDIX I. 


ELECTRICAL CONTROL AND CORRECTION OF BETATRON FLUX. 

The effect of applying D.C. bias to betatron-starting bars has already 
been considered in the main paper, and it is shown that considerable 
improvement in performance may be obtained. Other methods of counter- 

acting non-linearity in the betatron bars by auxiliary electrical circuits are 

available, and thetechniques involved are being developed at Metropolitan- 
Vickers Electrical Co., Ltd., in association with the 375 MeV. synchrotron 
magnet design. However, a brief review may be appropriate to this 
“paper. , 

One possible general technique is suggested by Kerst’s (1945) flux 
forcing circuit, which consists of an auxiliary winding around main poles 
and betatron bars, automatically ensuring the correct relation between 
annular and central fluxes. In application to a synchrotron the circuit 
must be ‘‘ switched off ”’ after betatron-synchrotron transition has occurred, 
to avoid excessive currents. Mechanical or electronic switching are 
possible solutions, if care is taken over transient effects. A more elegant 
modification is an electronic circuit in which, for example, any unbalanced 
voltage from a Kerst type winding is fed into an amplifier which in turn 
feeds a correcting winding on the betatron bars. Switching should then 
be relatively simple. 

Electrical control may also be used to give fine adjustment of the 
‘equilibrium orbit radius, an essential operational facility. Mechanical 
adjustment of the betatron air gap alone may be neither convenient nor 
sufficiently sensitive. (See for example, the air gap tolerance figures of 
‘§2(iii)). Anauxiliary bar winding connected across a variable inductance 
‘conveniently achieves the desired result. 


APPENDIX II 


375 MEV. SYNCHROTRON PERFORMANCE CHARACTERISTICS. 

In the main paper the design methods discussed were illustrated at 
various points by calculations based on the 375 MeV. synchrotron, and 
in particular it was determined (§§3 (ii) and 5 (i)) that a betatron bar of 
100 cm. effective length and 4200 sq. em. cross section, of the specified 
material and used without D.C. bias, would make at least 65 per cent 
betatron-synchrotron transition possible. In the specification finally 
adopted for the machine, the effective flux bar dimensions are 100 cm. 
length and 3300 sq. cm. cross section, and using this specification it was 
‘considered worth while to calculate the expected performance in some 
detail. Some of the more instructive results of the calculations are 
presented below. Since the machine will initially be operated at 300 MeV., 
some results for this energy are included. 

Using a hysteresis loop corrected for skin effect by the method of §6, 
the instantaneous circle motion of an electron, initially injected on to the 
‘equilibrium orbit, was calculated for various flux bar D.C. bias values, by 
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means of equation (26). The difference in skin effect for 300 and 375 MeV. 
operation is too small to affect the results appreciably. Fig. 5 shows the 
calculated orbits, and a superposed curve shows the synchronous orbit 
radius expanding towards 1, as the electron velocity approaches that of 


Electron Kinetic Energy (MeV) 
1@) Z a S 8 10 6 12, 14 16 18 


Orbit Radius (cms) 


~O 100 200 300 400 S00 
Magnetic Field at 12§ cms. Radius Gauss) 


Variation of orbit radii. 
375 MeV. synchrotron, 


a) Synchronous orbit. ; 
Bete electron orbit. Zero betatron bar flux bias. __ 
(c) Betatron electron orbit. 2000 gauss betatron bar flux bias. 


light. Betatron-synchrotron transition takes place at the intersection of 
this curve with a betatron orbit curve, and it would appear that transition 
conditions may be materially improved by flux biasing or other electrical 
correction methods. Further adjustment may be acheived by alteration 


‘ E 
SER. 7, VOL. 41, NO. 312.—JAN. 1950 


50 J.J. Wilkins on the 


of the resonator frequency, making the synchronous radius tend to some 
value other than 125 cm. 

The betatron output energy of the machine may also be deduced from 
tig. 5; e.g. with no flux biasing, electrons would strike the internal gun 
at about 12 MeV. energy. (An external target and internal gun are to be 
used in this machine.) 
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Volts per revolution diagram. * 
375 MeV. synchrotron. 
(a) Incoherent radiation loss. 
(b), (c) Total volts per revolution required for acceleration to (b) 300 MeV., 
(c) 375 MeV. 
(d), (e) Volts per revolution supplied by betatron effect, (d) 300 MeV., (c) 375 MeV. 
(f), (g) Volts per revolution derived from R.F. cavity, ( f) 300 MeV., (g) 375 MeV. 


To maintain a constant orbit radius the electron energy must increase 
by (27r¢/108)0H,/dt electron volts per revolution. Incoherent radiation 
of energy by the electrons (MacMillan 1945) is equivalent to an ener By loss 
(also in eV. per revolution) given by the expression (6/1077, )(E,/mgc?)4B3, 
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_ where E,is the total electron energy and myc? itsrest energy. Some coherent 
Saaintior loss also occurs, but assuming that the total charge associated 
with the electron beam will be about 1/100 (or less) of the theoretical 
maximum value 0-05 microcoulomb (Blewett 1945), this loss is likely to 
be comparatively small. Thus betatron effect and the R.F. voltage must 
supply to the electrons an energy represented by the sum of the two above 
‘expressions. 

Fig. 6 is an “eV. per revolution diagram ’’ showing the relative 
magnitudes of the above effects. The R.F. volts per turn (v) rises 
relatively slowly (adiabatically) to the maximum value vy. Synchrotron 
phase oscillation damping will considerably reduce electron loss during 
this period, giving more efficient overall electron trapping than that 
indicated by Goward’s (1949) methods. It is thus of interest to attempt 
an approximate estimate of the overall efficiency, taking some account 
of phase damping, but continuing to assume a simple rectangular R.F. 
voltage wave form. 

Bohm and Foldy (1946) and several other authors have discussed the 
phase oscillations of electrons in a synchrotron about the stable R.F. 
phase %, where ¢# is given by 


ih) Vee ee AY. a av. 4 (LE) 


When ¢ is zero, stable oscillations are possible up to a full 180° amplitude, 
but as & increases from zero, the permissible amplitude decreases and, 
moreover, the oscillations become assymetric about w#. Thus, when 
ab is 30° the maximum stable oscillation is from —40° to +150°, which is 
—70° to +120° about &% In order to make the following trapping . 
efficiency calculations pessimistic rather than optimistic, we consider only 
the smaller phase amplitude (such as 70° when ¢ is 30°), and denote this 
by 6,,- At all times during synchrotron acceleration the electron bunch 
must contain a range of smaller phase amplitudes from zero up to 8, 
where 0<6,,. 

At transition energy, v and % are approximately zero and 6,, is 180°. 
If some loss to the walls occurs at transition, then from the definition 
of f, in equation (14), the remaining electrons have phase amplitudes up 
to only (180° f,). After some synchrotron acceleration has taken place, 
and # is no longer zero, all electrons which tend to have phase amplitudes 
greater than @,, are considered lost. Thus, neglecting damping for the 
moment, of the electrons remaining in the beam immediately after 
transition a fraction wu will now remain, where (if @,,, is expressed in degrees) 


EW USO, cml meme oe ks cA (2) 
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‘This is essentially the calculation used by Goward (1949). 

It has been shown (e. g. Bohm and Foldy 1946) that synchrotron phase 
oscillations of small amplitude are damped in proportion to (KH, cos byt. 
‘The phase oscillation equation can only be solved for all amplitudes when 
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is zero (e. g. Twiss and Frank 1949), and, from the form of the resulting 
equation, the rate of damping appears to increase with the amplitude, 
approximating to E} over a wide amplitude range. Thus if the small 
amplitude damping law is assumed to hold for all amplitudes and % values, 
the tendency again should be to a pessimistic rather than optimistic 
result. Equation (2) is accordingly modified to the form. 


i (9/1 80f,)(E, cos /E,)*. sores vis hay eee (3) 
Fig. 7 
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Overall electron trapping efficiency. 
375 MeV. synchrotron. 
(a) Transition at 8-0 MeV, total electron energy. 
(6) Transition at 6-5 MeV. total electron energy. 
(c) Transition at 5-0 MeV. total electron energy. 
(At 0-75 em. off the donut centre). 
(d) Efficiency factor fs (neglecting phase oscillation damping). 


Fig. 7 shows overall trapping efficiency curves, determined by the 
data of fig. 6 and equation (3). Curves (a) and (b) assume transition 
exactly at the donut centre, while curve (c) represents the uncorrected 
transition of fig. 5, for which the betatron beam is displaced by 0-75 em., 
from the donut centre. For convenience of comparison curve (d) has 
been added, representing the efficiency factor fs of fig. 2 and showing the 
effect of neglecting electron oscillation damping. 
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As in Goward’s calculations, the rising part of an efficiency curve 
represents conditions under which loss of electrons during the synchrotron - 
phase is predominant, and the falling part is determined by the loss at 
transition. It appears that an R.F. voltage of about 5kV. peak will 
ensure 70 per cent overall trapping efficiency, without flux correction, 
and higher efficiencies are possible under adjusted conditions. Moreover, 
the effect of finite R.F. switch on time, which has been neglected, will 
be to give additional efficiency. 


IV. On the Interaction of Colloidal Particles.—II1. Some Theorems 
Concerning the Free Energy of Double Layers. 


By 8. Levine, Ph.D., 
Birkbeck College Research Laboratory*. 


[Received July 26, 1949.] 


SUMMARY. 


Making use of the Poisson-Boltzmann equation, a number of theorems 
relating to the free energy of the electric double layers of colloidal particles 
in dilute sols are developed. A general proof is given for the equivalence 
of two methods for determining this free energy, viz. by (i) using the 
fictitious process of charging the excess ions in the double layers subject 
‘to the condition of constant surface density of ions and (ii) imagining 
that the particles are charged by bringing ions from the interior of the 
dispersion medium to the surface. In the case of two infinite parallel 
plates, the difficulty that arises due to the infinity of the surface potential, 
when all the ions are supposed completely discharged, is satisfactorily 
‘explained. 

The present work is an extension of some of the results obtained by 
Verwey and Overbeek in their recent book on the theory of the stability 
of lyophobic colloids. Our theorem of the equivalence of the two methods 
of evaluating the free energy is more general than a similar theorem proved 
by Casimir in this book. We are able to derive the Casimir theorem 
by introducing special conditions. 


§1. INTRODUCTION. 


In their recent book, Verwey and Overbeek (1948) develop two methods 
of evaluating the free energy associated with the double layers of charged 
colloidal particles, namely by (i) charging the excess ions in the diffuse 
part of the double layers under the condition of constant surface potential 
and (ii) using a generalized Lippmann equation. A general proof of the 
equivalence of these two methods is given in the book by Casimir (1948). 


* Communicated by the Author. 
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In this paper we derive a corresponding theorem, where the distribution. 
of surface ions, and not the surface potential, is kept fixed, by means of 
an argument similar to that used by Casimir. Our equivalence theorem. 
can be derived without the use of the Poisson-Boltzmann equation 
(Levine 1949b) but the proof given here does depend on this equation 
(or rather a generalized version of it) and, in addition, only applies to: 
dilute sols. Our theorem is more general than that proved by Casimir, 
since it is valid if the surface distribution of ions is arbitrary (but unchanged 
during the charging process) and, furthermore, Casimir’s theorem can be 
deduced by assuming that the surface potential be constant during 
charging. 

Verwey and Overbeek apparently did not investigate our method of 
charging the ions in any detail because of the anomaly of an infinite 
potential (for the discharged state) in the case of two infinite parallel 
plates. This difficulty was also encountered by the author and, in fact, 
lead to an erroneous conclusion (Levine 1946) concerning the nature of 
the forces between the plates in dilute sols (an attraction was obtained 
at large separations, whereas the correct result is a repulsion). However, 
in §3, we are able to overcome this obstacle in a satisfactory manner. 


§2. ENERGY OF DoUBLE LAYER OF SINGLE PARTICLE. 

For the sake of simplicity we begin with a single colloidal particle of 
any shape immersed in a large volume of electrolyte. Let the surface 
charge arise from the adsorption of ions of type 1, charge e,: v is the density 
of such ions at position s on the surface. If ” is the total number of 


adsorbed ions, then 7= | sv, integrating over the surface (area) S of the 


particle. We imagine that all the ions in the system, including the 
adsorbed ones, are subject to the familiar Debye-Huckel charging process. 
Assuming that all the ions are being charged at the same rate, the fraction 
A will denote the stage in this fictitious process. Our problem is to: 
determine the free energy of the double layer at any stage A’<1. Now 
in an earlier paper (Levine 1946), we have already remarked that our 
method of determining this free energy is based on keeping the surface 
distribution of ions fixed when increasing A from 0 to A’. Thus v does: 
not depend on any intermediate value A<)’. It is also assumed that the- 
distribution v can be chosen arbitrarily, 7.e. we treat it as a thermodynamic 
variable. There is really no loss in generality if we assume that v has 
the form 


v=nv(X', 8), f ,v0(4, s) dS=1, { ccaipbrvehta Tae 


where 1(A’,s) is subject to the normalizing condition, but it is otherwise: 
an arbitrary function of A’ (which represents here the completion of the 
charging process) and of the surface position s. In this way n takes on 
the simple réle of a “* magnitude ”’ factor. 

The density of ions of species i(i=1, 2, ...., 8) in the bulk of the 
solution is assumed constant, independent of n, and will be denoted by no. 
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The average potential at position r in the diffuse layer of the particle 
is %, which is taken to be the same for the different types of ions. Casimir 
(1948) writes the Poisson equation for yb in the form 


V2 = —4rrp/D=Af (db), ee tsa ne eee Fo (2) 


where p is the average charge density in the diffuse layer, D is the dielectric 
constant of the dispersion medium and / is an arbitrary function. This 
becomes the Poisson-Boltzmann equation if we put 


fA‘s)= -->3 oy nie, exp (—Aegb/kT), Save wel s) 


where e; is the (full) charge of an 7 ion, k is Boltzmann’s constant and 
T is the temperature. A different form for f(A‘) has been chosen by 
Grimley and Mott (1947) when they take into consideration the hydration 
shell about each ion. The general conclusions in this paper do not 
depend on the particular form of the function f(A‘). The average 
electrostatic energy associated with the excess charge in the double layer 
at stage A is 


BA, nm; V)=Hez{ pr dSt+h[ ppdv. . . . . (4) 


The second integral in (4) is taken over the volume V of the outer diffuse 
layer, which may be assumed to extend to infinity, where the potential 
7% is zero. Then the corresponding electrical free energy at stage 2’ is. 
given by (see Levine 1948) 


i 
n)=2 | AHO Gc e WIL ON Ops mm 
ah 


The quantities sand p depend on QA, X’, v (and hence ”) and the coordinates 
r or s but to simplify the notation, we shall not indicate this explicitly. 
It will be clear from the text whether we are referring to the potential 
at position r in the outer layer or at position s on the surface. 

It is convenient to evaluate the quantity 


OF (A’, 2) as Fa OK(A, 2; A’) 
on oA on 
Since we keep v and hence w fixed when integrating with respect to A, 
n and ) are considered independent of one another, and we can differentiate 
under the integral sign with respect to n. The partial derivatives of the 
quantities % and v with respect to A or n will be denoted by a subscript; 
we note that v,=0. Substituting (2) for p into (4), we readily obtain 


B(A, m; XN’) ic 
See ray Va ex (b,v-+p,) ) dS— =f (b, f+ Ap, f) a ‘ (7) 


where f” if the derivative of f with respect to the argument A}. Now 
from (2) we have the relations 


Vb =frab+rapyf, — Ven=M fh + + + + (8) 


dx. oe ee (6) 
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Then use of (2) and (8) and application of Green’s theorem yields the 
following formule for the volume integral on the right of (7) 


DA DA / Ob, OW; 
8a il: (sb, V4, —b,, Vp) du=— os | é (¥, On, hm 7) ds. . (9) 
Here n, denotes the outward normal at the particle surface into the region 
of integration, which is the inward normal to the surface. 
Suppose that the dielectric constant of the colloidal particle is Di 
and y' represents the potential function inside the particle*. Since 
there is no charge inside the particle 


Vy =0, V2st—=0, V2 00 i ural ae 
and on the surface we have the boundary conditions 
0 Cyst 
=yl Dee De er 
p=, D ane D aa 4nvie, (11) 


We proceed to evaluate the various contributions to (6), which come 
from the surface integrals in (9). Interchanging the order of integration 
with respect to A and S and then integrating by parts, 


es Dd . Ob Seek D Of, Dee Ob na 
im), 7) bat = al [Yat |, Sta), Of epee 
(12) 


It is seen from (11) that 0/0n,=0 when A=0, although it is possible that 
#A0. On applying the boundary conditions (11), and noting that v,=0, 
the right-hand side of (12) can be expressed as 


D! Oust a 
Ne | Bae iss ne eee | | dS—e ad | 
1 _vaeaYn Aor I ip an, ee ey iF bv dS 


DI a : Oth, : 
+7), 2) ¥ One oe pate. Ge See ee a 


where ,_, is the (surface) potential at stage \’ and position s. This 
disposes of the first term in the surface integral in (9). The second term 
in this surface integral yields the contribution 


Dir , 2 (ee pieat: 
eae | bn ay (Se, ) =e | an 


8 


Di 1 Ot 
bv AS+ lk: ar{ WSs 


S n On, 
(13) 


to (6). If we now combine the three integrals in (127) and (13) which 
involve the function ~' and make use of the first of the boundary 
conditions (11), then a simple integration by parts yields the expression 


DI (# cyt at 
d T n I A 
rl Af (age Wag )dS=0, Bree Saget 


a “a 
In the case of a metallic particle, we put #'=constant. 
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‘This follows on application of Green’s theorem and of the equation (10). 
Also, we have two terms, one from (127) and the other from (13), which 
‘cancel out with the contribution from the surface integral on the right 
of (7). We finally derive that 

OBA oer ; 

eee === /\ ex] Panam dS. . : . . . (15) 


The formula (15) can be given a simple interpretation. When we 
bring up 6n ions, each of charge A’e,, to the surface, the increment in the 
‘density of ions at position s is 6y=v,5n—»,(A’, s) 8. The contribution 
to the change in electrical energy from the surface element dS is 
X'e,%,_.,,5v dS which is the electrical work done in bringing 8v ions to the 
element dS. In particular, if the density of ions is constant over the 
surface, then v9(A’, s)=1/S and (15) becomes 


oF, A, n) yy , 
oe ==) Cb 50(A 3 n), ey oe les, katte OSes (16) 


where y,,.(A’, 2) is the average potential on the surface. If we define the 
average potential in terms of the weight function v,=v,(\',s), then (16) 
‘can be interpreted as covering the general case (15). We now integrate 
with respect to » from 0 to n’, say, .and so obtain 


B,(N, 0’) =Ne, | RRC NGSGY | ON mtaeanres ine) ect LT) 
0 


where, in general, 
BavelA’s n)= Ip brawl’, s) ds. 2 : o . q (18) 


‘The relation (17), which was already implied in a paper by Derjaguin (1940), 
“was proposed without proof by the author in an earlier paper (Levine 
1948). The right-hand side of (17) is the usual formula from electro- 
statics for the work done.in bringing ’ ions, each of charge Xe, , by 
‘infinitesimal amounts from the solution to the surface of the particle, 
‘the potential of which is #,,,(A’, 2). 


§ 3. ENERGY OF TWO PARTICLES. 


We can readily extend the result (17) to a dilute sol containing P 
‘particles. Suppose that there are /, ions of type | on the surface of the 
pth particle. We imagine that these ions are transported from the solution 
to the surface at the same rate for all the particles and that the fraction 
€, which replaces the variable n, denotes the stage reached in this process. 
‘Thus n,=&n, will be the number of ions on particle p at stage €. The 
potential 7%, and the surface density of ions v, on particle p will be con- 
-sidered as functions of A, € and the position on the particle surface. The 
derivation of a simplified expression for the free energy follows the same 
-steps as for a single particle provided we replace e, by ,¢, for the pth 
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particle, differentiate with respect to € instead of m and sum over all che 
particles. The electrical free energy of the P double layers at stage A 
can be finally expressed in the form 
2 “1 Pty : 
Ne, 2 my | yare(r’, én’) dE=Ne, Z | pare(r’, m,) dn,» (19): 
p=1 0 p=1-0 
where #/2°°(A', n,,) is the average potential on particle p. 

In particular, we can apply this formula to two parallel plates, each 
of area A, at separation R. Since we are dealing with closed surfaces 
it is necessary to consider both sides of each plate and also the edges. 
The latter, however, can be neglected for sufficiently thin plates. In 
addition, we need not concern ourselves with the contribution to the free 
energy from the outer sides, since this will not vary with the separation R.. 
Making use of equation (9) in a previous paper (Levine 1948) the electrical 
free energy due to the excess charge in the two inner double layers can 
be written as 


BV, R)=2ANe, { O(’, », R) de, arcana oh 
0 


where A’e,D(A’, v, R) is the potential of the mean force acting on an 
adsorbed ion, due to the excess ions in the layers; the zero of this potential 
is taken in the interior of the solution. For finite plates, both v and 
@(X’,v, R) would not be constant over the inner faces because of the 
edge effects and we would have to introduce averages. However, this. 
correction is neglected here. The result (20) follows from the general 
theory based on the methods of statistical mechanics, which is developed 
elsewhere (Levine 1949b). If we now introduce the approximation of the 
Poisson-Boltzmann equation and identify ®(A’,v, R) with the potential 
function %,(A’,,),(p=1, 2), then the formula (20) is immediately 
derived from (19). In particular, if v is sufficiently small, @(X’, v, R). 
can be assumed proportional to v, in which case we have 


F(a’, ’, B)=Ad’ey’@(1',v', RB). 5... (208) 


This result was already suggested in an earlier paper (Levine 1948). 
When the linear Debye-Huckel equation* is applied to two (infinite), 
parallel plates, the potential between the plates is given by 


Arve, cosh Akyx on att 


A EY SH ae S 0,2 91), 
De, sinha’? “Dep < mes - 2) 


where # is measured from the median plane, so that we assume 
P(A, R/2)=@(A,v, R). We find that an additional term appears on 
the right-hand side of (207), namely —87Av’e?/Dx3R. This difficulty 
was met in a previous paper (Levine 1946) and a lengthy discussion by 
Overbeek (1946) and the author (Levine 1946 a) followed concerning this 


SETI" SE cTanreeeseer ee acini aaa Sis 
* We expand f(Ay) in powers of and retain the first term only, namely 


P(A) Ane. 
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discrepancy. Although the essential basis for the appearance of this. 
extraneous term was first proposed by Overbeek, no conclusive proof 
that (207) is actually the correct result was given. We are now in a 
position to show clearly the origin of this additional term. 

When all the ions are discharged, the formula (21) yields the 
condition 


Hb Ab (AS) = Save,/DatRaon,. say, bs...) (22) 
2=0 


throughout the region between the two plates and the fact that the 
potential (A, %)>0o as 1/A when A->0 creates a difficulty. Use of the 
Boltzmann distribution law yields the expression n? exp(—e,m/kT) for the 
limiting density (at A=0) of 7 ions at position w But if the short-range 
van der Waals forces near the particle walls are neglected, then 
necessarily the completed discharged ions should be uniformly distributed 
throughout the dispersion medium, which means that the condition 
lim Ais(A, x) =0 must be satisfied, in contradiction to (22). 


The analysis for a single plate developed in the previous section can 
be applied to each plate, provided the volume of the diffuse outer layer 


does not extend to infinity but to the median plane. If now lim Ap+40: 
A>0 


then the integral in (5) still converges but this is not true of each of the 
separate terms obtained by substituting (4) into (5). To overcome this 
difficulty of divergent integrals, we replace the lower limit of integration 
A=0 by A=e, where «¢ is a small positive number, and after integration 
let <>0. An additional term is obtained on the right of (12), namely * 


oe tira Ib yee (ee eam OL#) 2 e (23) 
dir 550 A>0 

which does not vanish, since 0/4/dn, behaves as A whereas 7 behaves as 

1/A when +0. The form on the right of (23) is obtained when we apply 

the condition that F(A, »)+0 as A>+0. We now need to introduce a 

corresponding correction to the right-hand side of our basic relation (17), 

which becomes 


B(N, 0’) =e ie ae n)dn—lim de, | ave(dsm) dn. . (24) 


The second term on the right of (24) was suggested by Overbeek (1946): 
although his explanation was not quite adequate. In the case of each 
of our two parallel plates, this new term becomes 


Maen lacs 


R 4 Av'*e? 
see sge as eon ols ee hea (2, 
=—Tim Ade, [4¥Q 5) a Dek (25) 


* We note that at the median plane, 0%/@n,=0, since the potential is a 
minimum. 
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which is precisely half the additional (extraneous) term appearing on the 
right of (207). The anomaly arises because we are considering two plates 
of infinite extent. This model cannot be used to describe two finite 
plates in the neighbourhood of A=0 since it prevents the diffusion of the 
excess charge of the diffuse layers between the plates into the bulk of the 


solution as A>0. The necessary condition lim A’=0 ensures that the 
A+0 


expression (25) always vanishes and hence our formula (17) applies to 
each of the finite plates. 

For two identical particles, symmetrically orientated with respect 
to the line joining their centres of mass, we shall denote the (weighted) 
average potential over each particle, which is defined in (18), by 
bayelA’, 2, R). It follows from (17) that the electrical free energy of the 
double layers of the two particles is given by 


F(X, 0’, R)=20'e, | WON a BY dae cl eae] 
0 
This is the simplified form that we set out to derive in this section. 


§4. Casrmmir’s THEOREM OF EQUIVALENCE. 


This theorem of the equivalence of the two methods, employed by 
Verwey and Overbeek (1948) to evaluate the free energy of the double 
layer, can be derived from our general results in §2. It will be sufficient 
to consider a single particle, since the argument is readily extended to an 
arbitrary number of particles. Also we shall first derive the Casimir 
theorem in the special circumstances when (16) is replaced by 


el) — Veit Weg ek ae 


On 
where (A’, n) denotes the constant potential at the surface, 7. e. it does 
not depend on the position s, but it is a function of A’ and». It should 
be noted that » is treated as an arbitrary parameter in (27). The following 
‘considerations show that (27) is valid for any ” in two particular cases. 

Suppose that we replace the general form (1) for v by 


v=n*(N/ (A, 8), [ r8O, 8) 0S lie 5 fs Poet Ieee 


such that in (18), %,_, is uniform over the particle surface. In these 
circumstances, it might appear that (27) applies, but we find that in 
general, this is not the case. For let the arbitrary function Vo(A’, Ss) in (18) 
be replaced by v¥(A’, s) but at the same time n be given any value. Then 
the surface potential #,_, will depend on » and, in general, also on the 
position s. If we now introduce n=n*(X’) as well as the distribution 
vj (A’,'s), then %,_, does become uniform over the surface. It follows 
that for this particular value of n (and distribution v&(A’, S)), Baye (A’, W*(A’)) 
can be identified with (A’, n*(A’)); in other words the relation (27) becomes 
valid if we susbstitute »=n*(A’). The problem is whether wy is 
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constant over the surface when the relative distribution is again v%(J’, s) 
but ” is arbitrary. One case where this is true is that of small n, when 
(2) reduces to the linear Debye-Huckel equation, with solution =n. 
where the function ¢ depends on the form of v*(A’, s) but not on n. For 
we can choose the distribution v%*(A’, s) such that the value of ¢ at the 
surface (at stage A’) is constant, independent of the position s. Denoting 
this constant value by ¢,_,, (27) will apply, with #(A’, n)=nd,_,.. 
However, if m is not small, the essential property that % be proportional 
to m no longer holds and (27) is not valid for a particle of general shape.. 
The other case arises when we are dealing with a spherical, cylindrical 
or plate-like particle, or with two parallel plates, provided the surface 
density of ions is uniform (equal to n/S) and any edge effects due to the 
finite size of the particles are neglected. 

Now in §2, we can treat \’ as a parameter, just like n, and we obtain 
the relations 


Vb NP by, Vb 


where the subscript A’ denotes partial differentiation. Using precisely 
the same method as in §2, we can derive quite generally 


af" ‘1 CK(A, n; 2’) 


mee ee drA=N'e, I. wb, vy AS=N'e,P(X',n), say. (28) 


It is observed that (28) immediately follows from (6) and (15) when 1 is 
replaced by X’._ When (27) applies, then #(A’, n)=y,_, and hence 


P(N, 2) =H, n) 5. | dS ]=0, Me oie (28:) 


since the surface integral in (287) is equal to n. 
Since it is simpler to treat the case where the potential at the surface: 
is constant (equal to (A’, n)), we begin with the relation 


A nv 
2 { Epis VY) d=Ne, | woumaging - bermeaphy: 
(0) 


which follows from (5) and (27). Remembering that A’ and n’ are con- 
sidered to be independent of one another, we differentiate both sides 
of (29) partially with respect to A’, make use of the relations (28) and 
(287), and then put A’=A. In this way we obtain 


5 B(., N= 65 [. w(A, n) dn rey I, (A, m) dn. . . (30) 
This is an identity in A and n’, and in particular, will be true a n’ is 
chosen that function of 1 which is defined by requiring that #(A, »’)=y5, 
a constant. This function is given by the inverse relation which we 
write as 

Wb att (As hp, peat a ~syeome. ow (BL) 
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If the number of surface ions is given by (31), the electrostatic energy 
.of the double layer can be expressed as 


B(A, (As Yo) s Y= Freon. PoE [ [Help = + (32) 


where the surface potential is constant (equal to %») when integrating 
over the volume V on the right-hand side. If we also substitute (31) for 
n’ into (30), then it follows that 


nA, Wo) 
|, [pph,,do=—e, | nA, #) )ay-+re, |” wild, m) dn. . (38) * 


‘The first term on the right of (33) has been written in place of 


NA, Wo) 
| (A, m) dn—eqihyn (ds Ho 


0 
‘Since ,=#(A, n(A, %o)), we have the relations 


fi.) dn= BT pasm dn—dons bo) | == 5 [n.d a 


0 
and hence (33) becomes 


5 | lel do—— sy [ dex [m(d, fy a |. . (34) 


Integrating with respect to A from 0 to A’, we finally derive 


[SI we. du=—X'e, {? n(X’, x) ss atti ey | n( (A, ws) dib. (35) 


For plate-like particles, the second term on the right of (35) vanishes 
and we obtain the formula derived by Casimir (1948). For a cylindrical 
particle, the solution of the linear Debye-Huckel equation yields the 
formulee 


(A, 2) = — pee n(A, who) = Drcottt ees (36) 
Kot Ky (Ang) ey K (AK a) 
where ” is the number of surface ions per unit length of particle, a is the 
radius of the circular cross-section and K, and K, are Bessel functions 
of purely imaginary argument. Since K,(Axga) and K,(A« a) behave as 
‘Ind and 1/Axoa respectively when A->0, the second term on the right of (35) 
again vanishes, although n—> oo (which is, of course, physically impossible). 
In the case of a spherical particle of radius a, use of the linear 
Debye-Huckel equation for % yields the formula 


p=poa exp[+Ac(r—a)]/r, 2... 1. (87) 


for the potential distribution in the doable layer, where r is the distance 
from the centre of the sphere. Thus 


Ne Da 
b(A, 2) = CMa at (A, fo) = aie (1+ Aayt) - + (38) 
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‘So that (A, Yo) 0 as A>+0. We now find that 


Po 
lim de, | nA, ¥) dp=Dap2/2, . . . . . (39) 
40° 0 
If we substitute (37) and the corresponding volume density of charge 


—(D/4m)V2p=—Drx2bi4a, . . . . . (40) 


which follows from (2) the Debye-Huckel approximation f (Ab) =AK2p, 
the left-hand side of (35) becomes 


—Da? an A exp (2AK_a) an| exp (—2Axkyr) dr= —(Darp2/2)A'kga. (41) 


However, it follows from. (38) that the first term on the right of (35) is 
SPADA LEA) Se Ghee, aot ee (42) 

No mention of this exception to the Casimir theorem is made in the 
book of Verwey and Overbeek. The omission by Casimir (1948) of the 
second term on the right of (35) is made at the end of his proof when he 
assumes (in his notation) dy(A=0)/dn=0, 1. e. the derivative of the potential 
inthe direction normal to the particle surface vanishes at A=0*. However, 
this condition is not satisfied for a spherical particle since it follows 
immediately from (37) that Hien (0b/Or),_g=—Y,/a. This derivative of 


the potential does vanish at the bree of a plate or cylinder. 

It is possible to obtain a generalization of the relation (35) for a particle 
of any shape, the surface of which need not be at a constant potential. 
On the right-hand side of (29) we replace the function #(A’, 7) by 
Bve(A', %) which is defined by (18), with the distribution »,(A’, s) chosen 
arbitrarily. Differentiating the modified form of (29) with respect to 
dN’, we again put A’=A. It is now assumed that »’ is chosen, such a function 
of A that the quantity ¢,,.(A, 2’) becomes equal to a constant yp, 
independent of A. This function will be denoted as 


Be ly Ae et ed ateasng Sebel yn (SL 1) 


which replaces (31). (Since #,,,(A, 2’) refers to stage A, it is defined by 
the weight function v,(A, s) in (18).) The operations performed on the 
right-hand side of (29) are precisely the same as those in the analysis 
above. Turning our attention to the left-hand side of (29), it is seen 
that we are charging the ions from 0 to A’, subject to a fixed surface 
distribution of ions, given by v=n'v,(A’,s). Also, an additional term 
now appears on the left of (30), namely re, ¥(A, »’), which will not vanish 
in the general case. When we reach the equation which is the 
generalization of (35), Beely, 


[Tf Welncend=—6 J P, gg (Pa tod) Ney | ngs) 
tim de, f'” Fae LEY Gel nt ean tga, ef (43) 
A—>0 


* Subject to fixed surface potential (45). 
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the integration with respect to \ on the left-hand side is being carried out. 
under the condition that, at each stage A, the average potential over the 
particle surface is fyye(A, Maye(As Yo)) =o, 2 constant. This result reduces 
to the Casimir theorem for a particle of general shape in the following 
special circumstances. The distribution v9(A, s) (at stage A) and the total 
number of surface ions ” are chosen in such a way that the surface 
potential is constant over the surface and equal to %. We can denote this 
distribution by v* (A, s) and the value of n by n*(A, yo), which now 
‘replaces (A, %)- Then it follows from (28) that Y(A, *(A, %o))=0 
and the relation (18) yields (at A’ =A) #,,.(A, n*(A, %o))=%;-,=%o. Hence- 
(43) takes on the same form as (35), which is the required result. 


§ 5. Discussion. 


We can suggest the following physical interpretation of our theorem. 
which is expressed in equations (5), (17) and (18). All methods of 
picturing the formation of the electric double layers of the particles. 
should Jead to the same value for the free energy provided such methods. 
describe thermodynamically reversible processes. Suppose now that the 
surface density of ions is treated as a thermodynamic variable, which we 
can change at will. Then the charging process introduced in the definition 
(5) of the free energy is carried out under equilibrium (reversible) conditions. 
and is therefore valid. Our second method which leads to the expression 
(17) for the free energy is also valid provided the electric work done in 
bringing the ions to the surface is against the potential of the mean force: 
acting on the ions. When we use the approximation of the Poisson- 
Boltzmann equation (the potential of the mean force being identified 
with the mean potential) and also restrict ourselves to very dilute sols, 
then this additional condition in the second method is apparently 
accounted for. 

It might appear that we have introduced some restriction by keeping 
the distribution (and hence number) of surface ions fixed when charging 
the double layers. However, this is not the case and the explanation 
depends on two theorems in (classical) statistical mechanics. Firstly,. 
it has been shown by Onsager (1933) that it is permissible to apply the. 
charging process to a sub-space of the positional part of the phase space 
of the ions. Such a sub-space is defined by a fixed (but arbitrary): 
distribution of surface ions but not, in general, by choosing the surface 
potential as constant. Secondly, the free energy has an extremely 
sharp minimum at the equilibrium distribution of surface ions. This 
latter property implies that the size of the particles is so large that we 
are dealing with a macrospocie system. By treating a large number 
of particles we have, in effect, macroscopic particle surfaces. (We ignore 
here the complication which arises when the adsorbed ions form part of 
the crystal lattice of the particles and which is briefly discussed elsewhere- 
(Levine 1949 a)). 
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The two methods of Verwey and Overbeek for evaluating the free energy 
lead to equation (35), with the second term on the right-hand side absent. 
In developing a physical argument to interpret the equivalence of their 
methods, these authors introduce the apparently arbitrary assumption 
that the charging process be carried out under the condition of constant 
surface potential. Previously (Levine 1946b) the author had raised the 
question as to whether the equilibrium conditions were being maintained 
in this charging method. In the case of spherical particles, for example, 
the number of surface ions n becomes infinite as A>0, when the surface 
potential is kept fixed. As a consequence, the two methods originally 
proposed by Verwey and Overbeek yield different results and the Casimir 
theorem of equivalence has to be modified by introducing the second term 
on the right of (35). 

It appears, therefore, that equilibrium conditions are not kept in the 
process of charging at constant surface potential. In the opinion of the 
author, little or no physical significnace can be attached to the artificial 
process visualized by Verwey and Overbeek. Nevertheless it is possible 
to obtain the correct expression for the free energy by this method, 
although not in all cases (e.g. spheres). The reason for the success of this 
approach is to be found in the nature of the theorems relating to the 
free energy, which are based on the use of the Poisson-Boltzmann 
equation and the restriction to very dilute sols. A number of these 
theorems have been developed in this paper, and others will be given 
later. Thus, we regard the equivalence theorem of Casimir as a convenient 
mathematical relation. 


In conclusion, the author wishes to thank Professor J. Th. G. Overbeek 
for helpful correspondence on a number of aspects in this paper. 
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ABSTRACT. 


The general problem of the detection of heavy ionizing particles, 
in the presence of y radiation and background noise from the amplifier, 
is considered. It is shown that the maximum possible signal to noise 
ratio is obtained if the input pulse is shaped to a form exp{— |t—T> |/7}.- 
The signal to noise ratio is also calculated for other common pulse shapes. 

A practical amplifier circuit is discussed which may be used to detect 
100 kv. energy protons in the presence of intense y radiation. 


INTRODUCTION. 


THE signal from an ionization chamber is a pulse of current which flows 
in an external circuit as the movement of ions induces charge on the 
electrodes. If free electrons are collected, and the counter geometry 
is properly designed, this time may be as low as 0-lys. which is less than 
the usual amplifier time constants. The total charge which passes is 
proportional to the energy of the initial ionizing particle. 

It is desired to amplify this current pulse and to distinguish it from 
background “‘ noise”’. This background noise is due to the superposition 
of many smaller pulses ; some due to electrons in the counter, others due 
to the spontaneous thermal agitation of electrons in the electrical circuit 
(Johnson noise) or the spontaneous movement of electrons in the 
amplifying valves (shot effect). 

The wanted current pulse, like all other current pulses, is transformed 
into a step pulse of potential on the grid of the first valve, by charging 
the stray capacity to earth of the ionization chamber and valve. The 
grid leak of this valve discharges the condenser after a long time, but 
the grid leak is usually made large in order not to produce appreciable 
resistor noise, so that this time is sufficiently long to be regarded as 
infinite. 

The amplifier will distort this pulse of potential. Part of this distortion 
may be deliberate, part caused by defects in the amplifier. It is convenient, 
however, to consider an ideal, non-distorting, linear amplifier and certain 
linear pulse-shaping networks which distort this and every other pulse. 
Then the problem is to decide upon the best way to distort this pulse 
in order to discriminate against the background noise. It is convenient 


* Communicated by the Author. 
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to “refer” the various sources of background noise to the input as an 
“equivalent noise charge’. This is a fictitious noise charge on the 
ionization chamber electrodes, which produces the same signal at the 
amplifier output as the noise. 

It will be shown that with a conventional resistance—capacitance 
pulse shaping network the best signal to noise ratio is obtained when 
the time constants of “ differentiation ” and “ integration ’’ are equal ; 
that an improvement of 20 per cent is possible by using a shorted delay 
time for shaping, and that a further improvement of 15 per cent is possible 
by using a network to shape the pulse to a form exp{— |t—T |/7}. The 
overall improvement of 35 per cent may be very important. 


ANALYTICAL EXPRESSION FOR SIGNAL TO NotsE RatTro. 


It will here be assumed that the duration of the current pulse is short 
compared with the pulse shaping time constants so that the input may 
be regarded as a function Qé8(t) of time—where 6(t) is Dirac’s 5 function 


+ 00 
defined by 65(¢)=0 for t40 and i o(¢) dt=1,. 


If this assumption is not true (as for example with ion collection in 
counters) the pulses from the ionization chamber will be attenuated 
more than the background pulses which is undesirable if a high signal to 
noise ratio is wanted. 

The potential charge on the first valve grid is then HW) where H(t) 


is Heaviside’s unit function. 


H(jje1 t>0, 
H(i)=0 t<0, 
and a = Oty. 


All pulses will be modified by the shaping network. If the input potential 
pulse is H(t), then the output pulse is a function f(t) of the shaping network. 
If the input potential pulse is 3(¢), the output pulse will be d/dt{f(t)}=f"(t)- 

The amplifier noise is due to the superposition of the effects of 
individual electrons in those amplifier elements which are at a low signal 
level—the input grid leak, the first valve and the first valve anode load. 
When the ionization chamber is used to count heavy ionizing particles in 
the presence of intense y radiation the secondary electron pulses in the 
chamber also cause a background. 

The noise is often expressed in terms of the amplifier bandwidth, 
here this practice is reversed, and the noise is expressed in terms of the 
output pulse shape, which is considered to be more fundamental in a 
pulse amplifier than the bandwidth. This involves a more fundamental 
consideration of the noise problem. 


EZ 
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The valve shot noise may be directly expressed as the superposition 
of pulses caused by individual electrons in the current stream. Each 
electron will produce a current pulse as it reaches the electrode: since 
the transit time (~10-® secs.) is small compared with the amplifier time 
constants, this may be regarded as a pulse of €6(t). 

Here can be distinguished two classes of noise. The valve anode has 
a resistive load R, so that a potential pulse Red(t) is developed across the 
load. The superposition of these pulses will be called class I noise. 
On the other hand, the valve grid has a capacitative load C, and a 
potential pulse eH(t)/C will be developed across the load. The super- 
position of these other pulses will be called class IT. noise. 

In practice it is desired to reduce the signal level at which accurate 
counting is possible. An arbitrary maximum rate of counting of noise 
peaks is set for each experiment, and the best amplifier conditions must 
be found in order to reduce the signal level at which this counting rate 
occurs. A useful concept for this problem is the root mean square noise 
fluctuation. If the noise current is I(t) the r.m.s. noise fluctuation o is 


given by y/(I(t)—I(#))?. 

If the rate of arrival of noise pulses is large compared with the 
frequencies passed by the amplifier, the rate of counting of noise peaks 
may be expressed in terms of o, the signal level V, and the amplifier 
bandwidth. (See equation (20) later). The problem of minimizing the 
signal level for a constant rate of counting is difficult, but the amplifier 
bandwidth only enters into the counting rate in a linear term, whereas 
o enters in an exponential term. A close approximation, increasingly 
accurate as V/o increases, will be to minimize o. A rough check shows 
that the pulses should be about 10 per cent wider than given by this 
criterion, but this would alter the minimum detectable signal by only 
1 per cent. ' 

The assumption that the rate of arrival of noise pulses will be greater 
than the frequencies passed by the amplifier will usually hold for shot 
noise from both the anode and the grid, and for thermal (resistor) noise. 
But the y ray background may not be sufficiently intense for the pulses 
to “pile up” in this way. In this case the y rays will be counted as 
separate pulses of height ke. 

Change of pulse shape will not then alter the rate of counting of 
background at all; so that the y ray background must be treated as 
distinct from all other forms of noise. If the pulses are small, their effect 
will be inappreciable, and if large, they will be similar to unwanted 
heavy ionizing particles. 

An intermediate case will often arise where a certain amount of pile-up 
occurs which. is difficult to treat theoretically, especially since the other 
forms of noise are added. 

If these limitations do not apply, we may proceed by finding the root 
mean square fluctuation from Campbell’s Theorem (Campbell 1909). 

Campbell’s theorem states that if a noise current I(t) is caused by the 
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superposition of v events per second, each causing an effect J(t) then the 
root mean square fluctuation o is given by 


=v | (fp dt. et ee a) 


The anode current develops a potential difference across a resistance 
and the pulse 8(¢) later is transformed into a pulse f(t) by the pulse shaping 
network. We therefore derive, following Schottky (1918), 


o2= Rey | P(t)? dt 
0 


=Ril,e | {PO} at 5 ss 2 = etme et 


Where o; is the r.m.s. anode current fluctuation. If we refer this to the 
lonization chamber input and put in the usual modification for space 
change Jimitation (North 1940), we obtain 


0-1 2ec? 1 8I 
=e ae ] aan 2 dt* 
Q Soe en + {f'(t) (t)} 


(where c is the total input capacity of valve and counter) . (3) 


for conventional valves with oxide-coated cathodes and a high amplifi- 
cation factor. 

The grid current produces a potential difference across a condenser 
and the fluctuation becomes 


ed=Ige{ fwd... . 2. 


The background secondary electrons in the counter will produce a 
background 


ok (ke) | {fly ae ee ee) 


Where Iy is the total ionization current due to y rays and ke is the 
average charge collected per electron. The noise in the input grid leak is 
(following Johnson 1925) 


o§=2kT/R | {fe} at era Po: <>" (6) 


which is small compared with the grid current noise if R is large, or more 
exactly if 
2kT/R <(I,+4l,)e. 
It can be shown that the anode resistor noise belongs to the same 
‘class’ as the anode current noise and is always small in comparison. 


* All r.m.s. fluctuations are now expressed in terms of a fluctuation of charge 
on the first valve grid. This is denoted by use of the suffix og. 
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Flicker noise cannot be classified in either class I. or class II., but is 
intermediate. The original process is frequency dependent and o* falls 
inversely as the first power of the frequency. It will only be appreciable 
for pulse widths greater than 1 millisecond. 

Induced grid noise is of the same class (class I.) as anode noise, but is 
always negligible in comparison. All the usual forms of noise can there- 
fore be represented in the form 


o§=Al (fp at+Bl {foprade . oh, aie 
0 0 
where 
A=(l;--hl,)e4 2U1/Ry =. ais = Sen 
0-12eC2 I, | 
Ba rt Bled) + ss + () 


: 0-12 eC? 
For a triode, B reduces to sae 
m 
OptrimMuM PULSE SHAPE. 


It is often of interest to obtain the maximum possible signal to noise 
ratio with no restrictions on the shape of the output pulse. This may 
be simply derived from the expression above by minimizing o§ for a 
constant peak value of f(t). f(t) must clearly be zero at t=0 and t= oo 
the first from initial conditions, and the second to make 


J. rep a 
finite. 
Moreover it should be a symmetrical pulse ; the integral { ce {f (t)}* dt 
0 


can be separated into 


Th co 
[rors] Oy ae 
0 Ty 


where f(t) is a maximum at t=T.' If f(t) is symmetrical the two integrals 
will be equal ; if f(t) is not symmetrical and 


TD ws 
| WP a> | {f@¥ at, 
0 aay 


then o will be made smaller by altering f(t) between 0 and T to make the 
pulse symmetrical. 


If f(t) is a function of (t/r) alone where 7 is some variable, then the best 


choice of 7 is given when the two integrals which add to produce o? are 
equal, 7. e. when 


e 


A] pare” (roy at. einen UA) 


img 
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if we can choose the form of the function f(t) then the best function is 
given by 

poO=exp.— i—Tal/ty s (9) 
and 72=B/A. 
This follows by minimizing 


[, Atrop act |” BUFO} diol]? 
with the auxiliary frac that ay 
| (= 
t=Ty f()= 


this is normalizing the noise to unit pulse height. Then using the 
Calculus of Variations we obtain 


Bd 
2 F(t 2? 
2f(t) ay {2f’(t)} (Euler’s formula), 


‘oni 


whence f(t)=exp{—(t—T>)/7} where for t>T) and f(t)=exp{(t—T>)/7} 
for t<T, by the symmetry property discussed above. There is a dis- 
continuity in f’(t) at t=Ty but no discontintinuity in {f’(¢)}. 

Clearly 


Af {fo}? ae—B | “(roy at 


in agreement with formula (8). 

From a physical point of view, it is of interest to see why this pulse 
shape is the best. The sharp pointed peak will not give use to excessive 
noise, since the extra height adds little to the area under the curve of 
{f()}?, whereas it adds eee to the signal. This pointed peak is 
the great advantage. 


COMPARISON OF PULSE SHAPES. 


It remains to be seen how much improvement can be obtained in 
practice by using this pulse shape in preference to other pulse shapes ; 
in order to do this it is necessary to consider the response of commonly 
used amplifier networks to a pulse H(t) or 6(¢) to obtain the shape f(t) 
and j’(t). 

The response of the amplifier to a function H(t) or 6(¢) may be readily 
found by the method of the Laplace Transformation. 


F(p)=[eef(e at Rene Bry eet (10) 


The transform of 6(t) is 1 and of H(é) is 1/p. 
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The system.function for the amplifier G(p) is written down by replacing 
the differential operator d/dt by p. Then the response to a pulse d(f) is 
the inverse transform of G(p) and the response to a pulse H(t) is the 
inverse transform of G(p)/p. The method is discussed by Elmore (1948) 
and Carslaw and Jaeger (1941). The method is here applied to obtain 
the signal to noise ratio for a variety of-pulse-shaping circuits. 

(1) Single differentiating circuit and single integrating circuit. 

The system function G(p)=72/(1+7,p)(1+72p) 


f()=(r9/7,)!-!™ fexp (t/7)—exp(t/7y) (721 72— 7) = As (11) 


when normalized to unity maximum amplitude (Elmore 1948). 
o2 is then a minimum for 71;=7T.=T=~/(B/A) and 


=3°7:4/(AB).< le sel ea 


rol S 


o=4/(AB), 


A similar response is obtained from using a critically damped tuned 
circuit with 7?=LC. 

(2) A shorted delay line in conjunction with a single integrating circuit 
gives optimum signal to noise ratio if the delay line time (for twice the 
line length) is equal to approximately 1-27 where 7 is the integrating 
time constant, and 

t= 4/(B/A). 


Then 
672-5 a (AB)o = bp 2 he see aed 


(3) A large number of integrating circuits limit the response of a multi- 
stage amplifier. After a large number of such circuits the response to 
a pulse 6(¢) tends to the shape 


f(t)=exp{—(t—Tp)?/7?} (Elmore 1948), . . . (14) 
where Tp=n7, and r=7,+/(n). 


Where x is the number of stages. When 1 is large, this is also the response 
to a function H(t) if this has been differentiated with a time constant 
small compared with 7. Pulse shaping with a shorted delay line and a 
large number of integrating circuits will also give approximately the same 
shape. Then 


ot==2-54/AB with 7=4/A/Bo.' 5... (15) 
The optimum case with 
J (t)=exp{-- [ta By | frp ee ee eet Le 
with t=1/A/B gives 
oa 247A Bai oe cowl ie seo cee), ALB) 
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These figures differ from those derived by Elmore (1948). Elmore: 
defines a pulse width 


8) 


Qa | (t—ty)2f(t) dt 
W= — 


0 
f(o) dt 
0 


where 


—+__ (17a) 
i TL) Ob 
0 


and derives the minimum noise—neglecting the grid noise (i. e. putting 
A=0) for a constant value of W. It is of interest to compute W for our 
various pulse shaping networks when the parameter f, has been chosen. 
to make o? a minimum. 


= 

I 

, 
el > 


For the optimum case 


For the Gaussian case W2— 


RS) 


For the RC—RC case 


= 
Il 
s 
-) 9 
We Bl le 


For the RC—Delay line case W?= 


Thus the Gaussian case has the advantage of giving a smaller pulse: 
width—according .to Elmore’s definition. This is in fact because there 
is less of a “ tail’ to the pulse. 

For most purposes it will not be necessary to produce a circuit to shape 
the pulse to the optimum. The improvement over a Gaussian or 
RC—Delay line case is only 15 per cent. For some experiments, however, 
this extra 15 per cent may inake all the difference. This arises when low 
energy protons are to be counted in the presence of intense y radiation. 
(Wilson, Collie and Halban 1949). A circuit has been devised (fig. (1)) 
which will shape the pulses to something which closely approaches the 
optimum shape. It is necessary, for such a circuit, that only linear 
networks be used; otherwise, although a pulse H(t) may be shaped to 
a pulse, f(t) a pulse 2H(t) will not be shaped to a pulse 2f(?). 

Briefly the operation of the circuit is as follows. The input pulse (1) , 
is split into two parts. One half is converted to a square pulse (2) and 
then into a nearly linear voltage rise (3). The other half is converted 
into an exponential rise (4). The first shape is controlled by the delay 
length and VRI, the second shape by VR2. The relative amplitudes 
are controlled by the feedback resistor VR3, so that the initial slopes of 
both waveforms (3) and (4)arethesame. The double triode V3 is connected 
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so that the difference between the two waveforms is taken as the output. 
The pulse shape after the peak is exactly the optimum shape, but that 
before the peak rises too quickly (theoretically it should begin to rise 
an infinite sime before the peak). With suitable choice of constants 


ota D1 A/ AB. see a tS em ee 


_ This computation made above of the optimum pulse shaping network 
assumes that the valve anode and grid noise are independent and constant. 
If the valve operating conditions are altered, the grid current noises 
can be reduced, at the expense of the anode current noise, by reducing 
the HT voltage (and hence the space current) and also by slightly reducing 
the filament voltage. The minimum noise is given when AB is a minimum 
(provided that the amplifier time constants are readjusted to the optimum 
values). If only valve noise and y background are important this means 
(AL,+1,)/Gm must be a minimum. This is a slight extension of a formula 
of Gillespie (1948). This assumes that the anode current noise is com- 
pletely space charge limited, and therefore only holds for a triode. 


Fig. 1. 


seit 


Pulse shaping circuit. 

If two different valves are to be compared, then the input capacity 
alters. A “‘ figure of merit’ for a valve in any given application is 
therefore 

ral Pll oteiL,).t abl inte eersteed Leen BE Lea 
Where C is the total input capacity of valve and counter (but not the 
feedback capacity due to Miller effect). 


CouNTING oF Norse Peaks. 
Before it is possible to compare these results with practice it is necessary 
to consider the rate of counting of noise peaks. 
The number of times that the noise waveform passes through a voltage 
V with positive slope is 
‘[exp{—V?/2V5}]{4 expected number of zeros of the noise waveform} 
(Rice 1944, 1945, Middleton 1948).. . 2... (20) 
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‘This is clearly the number of times the discriminator triggers. 
The expected number of zeros per second is 


| | ewona 
ih \ nr ra wiescarceeai Be Pee ee el) 
| J, Waar 


Where W(f) df is the noise power at frequencies between f and f-+df. 

For the cases considered this approximates to 0-6/7 where t= \/B/A. 
(The case of an integrating R-C circuit as the upper frequency cut off, 
-gives a divergent integral; this is discussed by Rice but will not arise 
in physically realizable circuits.) 

For most of our cases r= 10-8 secs. 

Then the number of zeros is 0-6 x 108/second. 

Then there will be less than 1 count per minute on noise peaks if V/V)>6, 
-t, €. if the signal is greater than 6 times the r.m.s. noise. 


PRACTICAL AMPLIFIER CIRCUIT. 


The above analysis has been carried out in order to design a counter- 
-amplifier system to count low energy protons in the presence of intense 
-y radiation for work on photodisintegration of deuterium (Wilson, Collie 
-and Halban, 1949). The valve chosen for first valve is an American 
6AK5 connected as a triode by strapping screen and anode. This reduces 
‘the anode current shot noise. The total input capacity of valve and counter 
is13 pF. The grid leak is made up only of the valve and insulator leakages 
-and is so high that no noise arises. 

The 6AK5 valves used were well aged, and then selected; the grid 
-current was then smaller than any other type tested ; the Mullard EF37, 
-although it has a low grid current, gives more anode current noise, and 
the EC91, or EF91, although they give little anode current noise have a 
higher input capacity and more grid current noise. 

Johnson and Johnson (1936) and Parsegian (1946) recommend pentodes 
for first valve to avoid the “ Miller effect” feedback. This feedback 
does not affect the signal to noise ratio, but it affects the gain and may 
result in accentuation of second valve noise. This can be obviated by 
a compensating positive feedback from the second valve anode to the 
first valve grid (Schultz 1946). We have, however, inserted an extra 
valve in the “‘ cascode ”’ connection (Hunt and Hickman 1939, Wallman, 
“MacNee and Gadsden 1948). The anode of the first valve feeds directly 
to the cathode of the second which represents a sufficiently low impedance 
(100 Q) to prevent appreciable amplification and feedback by the first 
valve. Yet the second valve noise is reduced by the high impedance 
-of the first valve in the cathode. This circuit also enables the first valve 
to be placed away from the subsequent valves—separated by a yard or 
:so of cable (fig. 2). 
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After amplification by a feedback ring a 100 @ line is taken to the main 
amplifier where the pulse shaping networks are situated. A shorted. 
delay line and resistance capacitance integration are normally used. it 
is possible, with selected valves, to make grid and anode current noise: 
equal with r=0-5 x 10-6 secs. with g,,=4 mA/V. Under these conditions, 
the equivalent number of ion pairs for the noise=o/c=300. 

Is is therefore possible to count pulses of ionization of 2000 ions pairs. 
above background, or 70 keV. 

When a y ray source is brought near the high pressure counter used, 
the value of 7 is reduced, to 0-3 10-8 sees. (below which loss of pulse- 
height might occur) and it is only possible to count down to about 90 keV. 
These values of the background, and the rate of counting of noise pulses. 
(equation (20)) are fully borne out experimentally, confirming that there 
are no other forms of background than the statistical ones considered. 
here. 

Fig. 2. 
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Amplifier input circuit. 
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ABSTRACT. 


This paper describes a method of finding the heights of occurrence of 
‘meteor trails from the radio echoes obtained by broadside reflection on 
a frequency of 60 Mc/s. The range of the reflecting point from the 
observing station is measured by the radio pulse technique and its 
‘elevation is determined by comparison of the signals received in two 
aerials mounted at different heights above a horizontal reflecting surface. 
Results obtained by this method during the Quadrantid shower on 1949 
-are in good agreement with those of visual observers. 


§1. INTRODUCTION. 


It is well known that a meteor, in its passage through the atmosphere, 
produces a column of ionization which can be detected by radio echo 
methods. Such methods have been developed independently by Hey and 
Stewart (1947) and Lovell (1948) in Great Britain; by McKinley and 
Millman (1949) Canada; and by Manning, Villard and Peterson (1949) 
United States. 

It has been shown by Hey and Stewart (1947), and by Lovell, Banwell 
and Clegg (1947), that radio waves of sufficiently high frequency (60 to 
75 Mc/s.) undergo specular reflection from the ionized column, the echo 
being returned from a small element of the trail situated close to the 
foot of the perpendicular from the station. A study of the durations 
of such transient echoes presents a number of interesting physical problems, 
and precise determinations of the heights of the reflecting points is necessary 
before progress can be made in their solution. 

Although the radio echo technique has been used previously to find 
the heights of occurrence of meteors in the atmosphere, the early methods 
depended on combined radio and visual observations (Prentice, Lovell and 
Banwell, 1947, Manning, Villard and Peterson, 1949), and were, moreover, 
not capable of high accuracy. A more recent method of improved 
accuracy described by Millman and McKinley (1949) can only be used 
on low frequencies (30 Mc/s.) where the reflection becomes non-specular. 

This paper describes a height finding system which operates on a 
considerably higher frequency, where the broadside reflection condition 


*Communicated by Dr. A. C. B. Lovell. 
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invariably holds. It is based on the spaced aerial method (Taylor and. 
Westcott, 1946), whereby the range and elevation of the reflecting point 
from a single observing station are determined by radio measurements. 
alone. The method was used to find the heights of 28 meteors during the 
Quadrantid shower of 1949, and has since been employed during the 
daytime showers of 1949 May to July. 


§2. PRINCIPLE OF THE METHOD. 


The apparatus operates on a frequency of 60 Mc/s. The range of the 
reflecting point of the meteor trail is measured by the radio pulse method, 
and its elevation is found by comparing the signals received in two aerials 
mounted at different heights above a plane, horizontal, perfectly reflecting: 
surface. 

For a scattering source of small dimensions at elevation % from the. 
station, the ratio, r, of the signals received in two such aerials is equal to 
the ratio of their effective power gains, G,(%) and G,(%), in this direction, 
and as these functions are calculable from the known parameters of the 
system, the relationship between % and r can be readily established. For 
a meteor trail, which subtends a finite angle at the station, the echoes. 
are due to the integrated effect of the radiation scattered from all parts. 
of the ionized column, and the evaluation of % in terms of 7 becomes 
more difficult. It can, however, be shown as follows, that provided the 
aerials are not too highly beamed, the amplitude of the echo received 
in either is proportional to its power gain in the direction of the reflecting — 
point. 

Consider first the case of scattering from a trail whose reflecting point 
lies close to the axis of a wide beam. The overall power gain G, of the 
combined transmitting and receiving aerial systems, can then be assumed 
to be constant over the first few Fresnel zones. It has been shown by 
Lovell and Clegg (1948) that under these circumstances the amplitude of 
the voltage fluctuations set up in the first stage of the receiver by scattering 
from a section of the trail lying symmetrically about the reflecting point 
and subtending an angle 26 at the station is 


+v +v 
V= KVR GH | cos 77/2v? dv+if sin 7/2? wh, ep ebiihy  § 


Where / is a constant depending on the parameters of the apparatus,, 
\ is the wavelength, R is the range of the reflecting point, « is the electron 
line density, and v=20,/(R/A). The amplitude of the signal received from 
the whole trail is obtained by taking the limiting values of the Fresnel 
integrals as v increases indefinitely. 

For accurate elevation finding it is necessary to make observations in 
directions where the gain of the spaced aerials varies rapidly with 0, so that 
toderivean expression for the amplitude of the signal received in either it is 
necessary to include G under the integral sign in (1). The value of V will! 
therefore depend in general on the radiation pattern of the aerial. However, 
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over a small angular range on either side of the reflecting point, G will. 
approximate closely to a linear function of 6, and may be written 


Gy+00G/00=G,+vdG/av. 


Where Go is the value of the gain in the direction of the reflecting point.. 
If the width of the beam is sufficiently great in comparison with the angle- 


subtended by the first few Fresnel zones this value may be substituted. 
in (1), which then becomes 


2 +0 +0 
\ =K(VR)*%a4 Go| cos (v?/2) dv+if sin (7v2/2) | 
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Plan of the apparatus (not to scale). 


The expressions under the integral signs in the second term of the right- 
hand side of this equation are odd functions of v, and the integrals vanish 
identically. The equation therefore reduces to (1), the amplitude of the 
echo is proportional to Gy, and the ratio of the signals received in the two 
aerials of a spaced system is equal to the ratio of their gains in this. 
direction. 


§3. APPARATUS. 


A plan of the apparatus is shown in fig. 1. The hut A contains a 
transmitter and receiver of conventional design. The transmitter works. 
on a frequency of 60 Mc/s. (A=5 metres), and radiates pulses of 8 y/sec 
duration and peak power 100 kw. The amplification of the receiver is 
linear, and its sensitivity is such that a signal of 2 .V. at the input doubles. 
the signal to noise ratio. The transmitting array, T, floodlights the 
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northern quadrisphere of the sky. It consists of a horizontal folded 
half-wave aerial at a height of $A above the ground, with a reflector 
adjusted to give minimum backward radiation. 

The two receiving aerials R, and Ry, are similar to the transmitting 
_array, and are mounted on posts 20 ft. apart at a distance of 200 ft. from 
the hut A. In front of these aerials is a plane macadam square, pqrs, 
with sides 100 ft. in length which can be covered with wire mesh to form 
.a perfectly reflecting surface. The heights of the receiving arrays above 
the square can be varied from 3 ft. to 18 ft., but they are normally 
maintained at $A and #A. The aerials are connected to the two stage 
preamplifiers P, and P, housed in the small hut B, and the signals pass 
from here through a common coaxial cable to the receiver. The 
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Theoretical relation between ratio r and elevation yx. 


preamplifiers are suppressed alternately at a frequency of 25 c/s. by 
multivibrator pulses supplied from the receiver. The same multivibrator 
pulses are applied to the a-plates of the A scope display tube, so that the 
echoes from the two aerials are displayed side by side. Their relative 
amplitudes are compared visually. 

For purposes of calibration a fourth aerial C is mounted at a height 
-of 20 ft. in the centre of the side of the square opposite to the receiving 
aerials. A duplicate display tube D with its associated time base and 
power circuits is provided in the hut B. Before each set of readings the 
aerials R, and R, are brought to the same height above the reflecting 
surface and a signal generator is fed into the aerial C. The gains of the 
preamplifiers are then adjusted until the signals received in R, and R, 
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are of the same amplitude, and the aerials are then returned to their: 
original heights of $A and 3A. . 

Fig. 2 shows the ratio r of the effective aerial gains plotted as a function: 
of the elevation, 4. As the aerials fulfil the conditions required by §2, 
this curve can be used to obtain the elevation of the reflecting point of 
a meteor trail from the observed value of 1. 


$4. LimrratTions AND Errors. 


(i) Error in Measurement of the Ratio r. 

With the apparatus in its present form the most serious inaccuracy 
arises from the uncertainty in the estimation of the ratio r. To find the 
magnitude of the error involved artificial echoes were displayed on the: 
cathode ray tube, and the apparent visual ratio was compared with that 
obtained by photographic measurement. The results showed that when 
the amplitude of the larger echo exceeds that of the noise by a factor 
of 2-5 or more, the ratio can always be measured to within | part in 10, 
and that for signals differing in amplitude by less than 20 per cent accuracy 
of 1 part in 20 is possible. For values of % greater than 20° the consequent 
error in the height is never greater than +4 per cent. 


(ii) Ambiguity. 

The vertical polar diagram of the higher aerial has two lobes with a 
zero at 41-8°, so that there are two values of % corresponding to each 
value of r. One of these usually gives an improbable value for the height 
H, and can be rejected. For a trail at elevations between 33° and 48°, 
however, there are invariably two possible values of 4%, and the method 
cannot be used. 


(iii) Error in Range Measurement. 
With the present apparatus ranges can be measured to within +2 km. 
This corresponds to an error in height measurement of +1 km. 


(iv) Imperfections in the Reflecting Surface. 

The method depends on the fact that the reflecting surface acts as a 
plane mirror with a reflecting coefficient of unity. It is well known that 
for this condition to obtain the surface must behave as a perfect plane 
reflector over an area comparable with that of the first Fresnel zone. At 
the lower limit of elevation (20°) the distance of the geometrical reflecting 
point from the higher aerial is less than 50 ft., but the first Fresnel zone 
extends outwards to a distance of 180 ft., and a considerable part of its 
area lies outside the surface pqrs in fig. 1. It is therefore important that 
the ground shall have good reflecting properties for some distance beyond 
the edges of the square. 

For 120 ft. beyond the sides pq, qr, and rs, the earth is at the same level 
as the prepared surface, and is plane to within A/16. The average values 
given by Burrows (1937) for the dielectric constant and the electrical 
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conductivity of soil are 10 and 5x10? e.s.u. respectively. This implies 
that for waves of frequencies of 60 Me/s. its refractive index, n, is equal 
to 10—i5/3. Substituting this value in Fresnel’s equation we find that 
for horizontally polarized waves the reflexion coefficient at a glancing 
angle x is given by 


sin p— 1/(10—7 5/3 cos? ) 
p™ sin b+ / (10--25/3 cos? p) ° 
In directions within 30° of true north the geometrical reflecting point 
lies outside the area pqrs only when ¢% is less than 7°, or when p is greater 
than 0-94. The modification in the vertical polar diagram due to this 
small departure from perfectly reflecting conditions is inappreciable at 
elevations greater than 20°. We may therefore assume that over the 
range of elevations with which we are concerned the curve in fig. 2 is 
unaffected by the limited area of the prepared surface. 


(v) Coupling between the two receiving aerial systems. 

To measure the degree of coupling between the two receiving aerials, 
a signal generator was fed into the aerial C (fig. 1), and the sigaa! from 
the aerial R, was permanently suppressed. The amptitude of the signal 
received in R, was measured by a microammeter connected across the 
output of the receiver. The aerial R,, which was initially at a 
height of 3 ft., was raised slowly to a height of 18 ft., and was then dis- 
connected. During this operation there was no perceptible change in the 
reading of the meter. The coupling between the two aerials must therefore 
be negligible. ‘ 


(vi) Change of aerial impedance uith height above the ground. 

The method of calibration described in §3 assumes that the impedance 
of the receiving aerials is independent of their height above the reflecting 
surface. To measure the variation in impedance with height a signal 
was radiated from the aerial C and received in the aerial R, in the way 
described in (v) above. It can readily be shown that if the impedance 
of R, is independent of its height, h, the strength of the signal received 
will be proportional to sin xh, where « is a constant. If, however, the 
mutual impedance between the aerial and its image is appreciable, the 
signal will show fluctuations in amplitude about the sine curve as h is 
varied. 

On varying h from 3 ft. to 18 ft. these fluctuations were observed, but 
their amplitude corresponded to less than 0-5 per cent of the total 


amplitude of the signal. The change in impedance with height is therefore 
trivial. 


§5. RESULTS OBTAINED DURING THE QUADRANTID SHOWER OF 1949. 


This method of height measurement was tested during the Quadrantid 
meteor shower of 1949. The periods of observation were from J anuary 3, 
03" 08" to 04740", January 3, 13"00™ to 13249", and J anuary 4, 
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03" 08™ to 04" 30™. At these times the ground was covered with a layer 
of wet snow to a depth of 5in. To find the effect of this snow on the 
reflecting properties of the surface, the experiment described in §4, 
section (vi) was repeated during the early morning of January 3. No 
deviation from perfectly reflecting conditions was detected. 

The mean echo rate during the periods of observation was 14 per hour. 
Previous experience with similar apparatus of comparable sensitivity 
has shown that the rate when no shower is active is two per hour. It 
may therefore be assumed that the majority of the echoes observed were 
due to meteors from the Quadrantid radiant. Out of a total of 41 echoes 
two were at elevations less than 25°, nine gave ambiguous heights, and 
in two cases there was uncertainty in the measurement of the ratio. 
The remaining 28 measurements gave heights between 80 km. and 107 km., 
with a probable error of less than 5 km. in each case. The height distri- 

- bution is shown in fig. 3. 


Fig. 3. 
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Height distribution for Quadrantid shower, 1949. 


$6. RESULTS OBTAINED DURING THE DAYTIME SHOWERS oF 1949, 
May To JULY. 


The method has recently been used to measure the heights of meteors 
of the summer daytime streams of 1949. The apparatus was slightly 
modified to give greater accuracy of range measurement, the pulse 
recurrence frequency being increased to 150 c/s., and the time base 
expanded to cover only 250 km. During the period 1949 May to July 
watches were maintained for a total period of 8-26 hours, and 239 echoes 
were observed. Of these 166 gave unambiguous heights. The height 
distribution is shown in fig. 4. 


$7. Discussion OF RESULTS. 
The results obtained by this method show that it is capable of giving 
G2 
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values for the heights of approximately 70 per cent of the total number 


of meteors observed. 
The beginning and end heights obtained visually by Prentice for trails 
of the Quadrantid meteors are shown in. Table I. 


TABLE I. 


Heights of Occurrence of Quadrantids obtained fromy 
visual measurement. 


Visual No. of Meteors Mean Beginning | Mean End 
Magnitude Observed. Height (H,) Height (H,) 


108 km. 84 km. 


0 

yy F 105 km. 87 km. 
4 103 km. 91 km. 
6 101 km. | 94 km. 


Mean Height of Occurrence eal =97-8 km.. 


Fig. 4. 
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Height distribution for daytime showers, 1949. 


It is evident that the heights shown in this table are confined within 
somewhat narrower limits than those obtained by the radio echo method 
but this is to be expected, since the values shown here represent the means 
of a relatively large number of observations. The value of the mean 
height given by Prentice is in good agreement with that obtained by the 
radio technique. 

Although it is not possible to obtain a direct comparison between the 
results obtained by this and other methods for meteors of the daytime 
streams, the height distribution shown in fig. 4 is of a form which would. 
be expected for meteors of medium velocity. 


Measurement of Heights of Reflecting Points of Meteor Trails 85 
§8. CONCLUSION. 


We may conclude, from the results listed above, that the method in its 
present form is capable of measuring the heights of the reflecting points 
of meteor trails to within +4km. The discussion in §4 indicates that 
by far the most serious source of error is the measurement of the ratio of 
the echoes. This may be eliminated by the addition of automatic photo- 
graphic recording, which will also reduce the range error. It is estimated 
that with the addition of this facility, a probable error of +1 km. will 
be attainable. 
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VII. On the Mass of Charged Particles of the Cosmic Radiation. 


By C. Franzinerti, D. in F., 
H. H. Wills Physical Laboratory, University of Bristol*. 


[Received October 31, 1949. ] 


[Plates I. & IL.] 


SUMMARY. 


The mass of charged particles, produced by cosmic radiation at an 
altitude of 11,000 ft., and brought to rest in photographic emulsions, 
has been determinedby ‘ magnetic-deflection ’ experiments. The method 
is applicable to particles, of any type, of which the lifetime is greater 
than 10-?° sec. 

The results show that at an altitude of 11,000 feet the particles are 
composed of different types in the following proportions :— 


ut-particles, 8-5=1-5 per cent; p—-particles, 8+ 2-4 per cent ; 
a+-particles, 4:0*1-0 per cent ; 7~-particles, 5-4=2-0 per cent ; 
Protons, 48+4 per cent ; deuterons, 16+2 per cent ; more massive 
positive particles, 11+2 per cent. 


The observations on the tracks of 360 particles give no evidence for 
the existence of other kinds of particles sufficiently stable to be detected 
by this method; and they show that if such particles exist, they are very 
rare. 


§ 1. INTRODUCTION. 


SOME time ago it was pointed out by Powell and Rosenblum (1948) that 
with a suitable arrangement of photographic plates placed in a magnetic 
field it is possible to measure both the momentum and the residual range 
of a charged particle, and thence to determine its mass. It was proposed 
that two plates should be placed parallel to one another and perpendicular 
to a magnetic field, the two emulsions being face to face. A charged 
particle leaving one of the emulsions at a small glancing angle, ¢, and 
with sufficient energy to reach the opposite emulsion, will be recorded 
in both plates. 

Owing to the magnetic field, the particle will move, in the space between 
the plates, along a spiral path lying on a circular cylinder of which the 
radius, p, will depend on the component of the momentum of the particle 


*Communicated by Professor C. F. Powell. 
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in the plane perpendicular to the field. From the angle between the 
directions of motion of the particle, at the points of exit and entry at the 
two emulsions, p can be determined. If the particle stops in the second 
emulsion, its residual range, 7, can also be measured. Inserting the 


appropriate values of p and r in the equations for the momentum, ‘p, 
and the energy E ; viz.: 


set eZ Hp 

Pagel ira) Saranp bovsid mere touts aee) whi 
and 
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the mass of the particle, m, can be determined if its charge, Ze, is known. 

It was suggested by Powell and Rosenblum that this method would 
give more precise values for the masses of the 7 and p-mesons than can 
be obtained either by means of grain-counts or by observations on the 
scattering of the particles. At the time, such improved measurements 
were of great interest because of their bearing on the nature of the neutral 
particle emitted during the decay of the 7- meson, but shortly afterwards 
the announcement was made of the artificial production of 7-mesons 
at Berkeley. It seemed clear that these experiments would give much 
more favourable conditions for the determination of the mass of the 
particles, and the cosmic-ray experiments were therefore temporarily 
abandoned. 

Interest in the method has been recently revived by reports, from 
several laboratories, of observations which have been taken to indicate 
the existence of other types of mesons in addition to the z- and p-particles. 
Further, some experimenters (Alichanian et al. 1948) have suggested 
that these other types occur with a frequency comparable with that of 
the z- and p-mesons. The question of the existence or non-existence 
of such particles is of such decisive importance for the development of 
nuclear physics. It was therefore decided to proceed with the magnetic 
deflection experiments, both because they appeared to be capable of 
greater precision than other methods employing photographic plates, 
and because they have the advantage that any secondary effects associated 
with the arrest of the particles in the emulsion can be observed. The 
second of these features is of particular importance because it is reasonable 
to assume that any other types of mesons must, like the 7- and y-particles, 
be unstable ; and any positive particles should therefore give distinctive 
evidence of spontaneous decay. 

A further important advantage of the method is that the time of flight 
of a particle in the apparatus is only of the order of 10~1° sec., a period 
much shorter than that necessary in experiments employing Wilson 
chambers and counters. It is therefore possible, in principle, to determine 
the mass of particles, by the photographic method, even although their 
existence is so transient that they would decay “in flight ” in traversing 
apparatus employing other methods of detection. 
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The present paper describes an application of the method to the 
determination of the mass of charged particles observed in photographic 
plates exposed at the Jungfraujoch High Altitude Research Station. 


§ 2. DerarLs oF EXPERIMENTAL PROCEDURE. 


Apparatus. The experimental arrangement employed is shown 
schematically in fig. 1(a). The plates were carried by a suitable plate- 
holder, (B), which was placed between the pole-pieces of an electro-magnet. 
A metal housing, (A), was fitted around the pole-pieces in order to maintain 
the plate-holder in its correct position during the exposure, and to make 
it possible to evacuate the space between the plates if it was found to 
be desirable. 

The plates employed were coated with Kodak NT4, or Ilford G5 
“* nuclear research ’’ emulsions, 504, 100 or 200 thick, and with linear 
dimensions 34314 inches. Emulsions thicker than 2004 have not been 
found convenient for the investigation, because of the existence of 
distortion in the processed plates. The most important advantage of 
using thicker emulsions is that there is an increase in the number of 
tracks which stop in the emulsion, and which are therefore measurable. 
The results given in a later paragraph show, however, that the trajectory 
of a proton with a residual range of 1000u has a radius of curvature 
of ~20cm., in the conditions of the experiment. The corresponding 
value of the deflection, if the path in the air-gap is ~2 cm., is only 6°, 
and such a small deflection does not allow a precise measurement of the 
mass of the particle. The gain in the number of tracks obtained in a 
given exposure with thick emulsion would therefore be balanced by a 
reduction in the average value of the precision of the determinations 
of mass. These considerations do not apply if the observations are 
confined to the particles of small mass, such as 7- or »-mesons, which are 
more deflected. It was important, however, for the purposes of the 
present investigation, to make measurements on particles of all types. 

Any small movement of the plates during exposure would have made 
it very difficult to correlate the two tracks produced by a single particle 
and would have greatly reduced the accuracy of the measurements. 
The plates were therefore attached to the holder with ‘ Kodaflat ”’ 
liquid rubber. The rubber was poured on the frame of the holder in 
sufficient quantity to make a layer 300, thick, and was allowed to dry 
for 12 hours. The plates were then pressed against the adhesive surface, 
and all possibility of their relative movement was excluded. 

Magnetic Field. A field of 29,180 gauss, ina volume of z X (3-5)? x 0-7. em.3, 
was maintained byan electro-magnet, the necessary power being 600 watts, 
at a current of 8-5 amps. The field was constant during the exposure 
to within 300 gauss, apart from ‘edge effects”, which were only 
appreciable in the region extending 4mm. from the edges of the pole- 
pieces. The variation of the field with the current is shown in Table I., 
and it will be seen that it is very small for currents between 8 and 9 amps. 
During the exposure, the changes in current never exceeded 0-25 amps. 
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TABLE I. 
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Variation of field-strength with magnetizing current. 


I (amps.) 8-0 8-2 8-4 8-5 8-6 8:8 9-0 
H (gauss) 28,900 29,020 29,140 29,200 29,250 29,340 29,440 
Bigs 1. 
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water-cooling 
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(3) window 


plate holder * 


Experimental arrangements : (a) during the exposure, 
(b) X-ray marking device. 


The field was determined by comparison with a known field produced 
by a permanent magnet, kindly made available by Dr. S. Rosenblum and 
Mr. Tsai of Bellevue Laboratory in Paris, using a fluxmeter of high 
precision. A cooling system maintained the plates below 20° C. 


Correlation of Pairs of Tracks. A precise knowledge of the geometrical 
relationship between the two plates during exposure was required in 
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order to correlate the pairs of tracks produced by each particle, and this 
was obtained by the following method :— 

The holder containing the plates was exposed to an X-ray beam, 
passing perpendicular to the plane of the emulsion; (see fig. 1(6)). 
Between the plates and the X-ray source was a screen of lead a few mm. 
thick, and provided with a narrow slit 8 cm. long and a few microns wide. 
With this arrangement it was possible to “ print’ a grid of lines on both 
the plates, as indicated in fig. 2, in such a way that corresponding points 
on the grids represented the same x, y, coordinates. The edges of the 
slit consisted of two cylinders of tungsten, 2 mm. in diameter and it was. 
made sufficiently narrow by pressing a piece of thin paper between the 
two cylinders during assembly. The black lines produced in the processed 
plate by X-ray exposure had a width of 20, or even less. 


The rectangular grid of lines on the two plates greatly facilitated the 
correlation of the tracks. The two plates of each pair were first scanned, 
and maps were constructed showing the track of all particles crossing the 
surface of an emulsion, and of those which stopped in one of them. From 
the position, direction, and angle of dip, ¢, of a track, T,, which stopped 
in one emulsion, (see fig. 3) it was possible to determine a region, R,, 
in the opposite plate, in which the continuation of the track should be 
found. If the density of the tracks in P, is low enough, only the 
true continuation of T, will have the correct ionization, angle of dip, 
and direction, consistent with it being the partner of track T,. The 
probability that the track produced by another particle will, by chance, 
satisfy the same requirements is discussed in §4, and is shown to be: 
very small, 
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Processing of the Emulsions. Great care was taken in processing the 
plates in order to reduce the distortions of the emulsion to a minimum.. 
The following procedures have been found satisfactory : 
eee RePee Ne PREIS CODA eo ieet Heth ENO rin (Pn ett. .  t 
Time 
100n 200u 
oo ee a? Wo Ee a ee Pa ee 
For Kodak NT4. 


Chemicals Temperature 50 uw 


Ts Co ae 0) \ ¥8> Ge 30m. 40m. 60m.. 
Hardening Stopping Bath 5 10m. 15m. 20m. 
Fixer (Hypo 30 per cent Solution 

in weight) a5 3/2 of clearing time. 
Washing in water 5 3h. 4h. 5h, 


For Ilford G5. 


Developer 3° C2 —— — 35m. 
16 ce. Azol Dev. 

10 ce. 1 per cent K Br. Sol. LspC; — — 200m.. 
800 cc. H,O. J 

ee BECP a ee \ As for Kodak NT4. 


§ 3. MeTHOD OF CALCULATION OF Mass. 

According to Camerini and Lattes (unpublished), the range-energy 
relationship for protons, in a photographic emulsion, can be represented 
by a simple power law for values of the range greater than 200, : 

EH=KR’"*, 
where E is the energy of a proton of range, R. Further, the ranges: 
R(v) and r(v) of particles having the same initial velocity, v, but with 
masses M and m, and charges Ze and ze, respectively, are related by the 
equation 
R(v)=r(v) . M/m(2/Z)P?—C, 

C being constant*. Since C is negligible compared with the ranges of 
particles considered in the present experiment, it may be neglected. 
It follows that the relation between the energy « and the range, r, of a 
particle of mass m, and charge ze, is given by 

NL Cy LW Be PAT oe Rete RA eee 
M being the mass of the proton, and ¢ the elementary electronic charge. 

The equation (3) is the explicit form of the equation (2). Solving the 
simultaneous equations (1) and (3), we obtain 


ma a(s fy jie A enn dig pest nhros,e (4) 


* (The constant C takes account of the capture and loss of electrons in the- 
last part of the range, a process which affects differently particles of different 
charge. See Bethe and Livingston (1937).) 
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on 
‘ 2 2—n 
where = ee and A= 13 min 2K : 
2—n CA a 
For ranges shorter than 200p, the range-energy relation given above 
is no longer valid, and the function can therefore only be evaluated 
numerically from the experimental points. If the logarithm of the 
residual range of a particle is plotted against the logarithm of its measured 
radius of curvature, it follows, from the equation (4), that the particles 
with equal masses should be represented by points lying on straight 
lines, of which the slope, #, is given by 


tan p=n/2. 


There will be a departure of the observations from this linear relationship, 
however, in the case of particles of short range. 


§4. Sources OF ERRorRs. 


The most important errors in the experiment are due to: (i) distortion 
-of gelatine ; (ii) scattering in the air-gap ; (iii) scattering in the surface 
layers of the emulsion ; and (iv) spurious coincidences. 

(i) Distortion of gelatine. 

Serious distortions of the emulsion extending over large areas, which 
can change the coordinates of a point by more than 10u, have been found 
to be confined to regions near the edges of the plate. The affected area 
extends inwards from the edges of the plate for 2mm. in the case of 
emulsions 504 thick ; and for 1 cm. for 200, emulsions. On the other 
hand, local distortions, which may considerably alter the direction of a 
track, sometimes occur if special care is not taken in processing. Local 
distortions are very often accompanied by differences in the general 
background of grains; and they may sometimes be detected, with the 
naked eye, in the form of manifest imperfections in the smooth surface 
of the emulsion. 

(ii) Scattering in the air-gap. 

In crossing the gap between the plates, the particle may be scattered 
by the gas filling it. The influence of such effects on the measurement 
of the curvature of the tracks of particles in magnetic deflection experi- 
ments, has been discussed by Bethe (1948). His calculations, when 
applied to our particular experimental conditions, lead to the result 


given in Table II., which shows the values of 22, the ratio of the radius 
Ps 
of curvature of a track due to the magnetic field, to the average apparent 


; 8s 
radius of curvature due to scattering. PE <a , where dg and 6, represent 
PS H 
the angular deviations dueto scattering and the magnetic field, respectively. 
(ii) Scattering in the surface layers of the emulsion. 
‘The determination of the direction of motion of a particle, at the exit 
from one emulsion and at entry into the other, is subject to errors as a 
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result both of scattering by the atoms of the emulsion, and of the dis- 
placement of the individual grains in the track from the true line of the 
trajectory, The second of these effects has been studied by Goldschmidt 
et al. (1948), in experiments with C2 Ilford plates, and has been found to 
be less than three per cent of the real scattering for tracks of range less than. 
2000. The error from “ spurious scattering ’’ can therefore be regarded 
as negligible. 

The errors due to scattering in the surface layers in the emulsion 
depend on the method of defining the direction of the track. If we choose. 
to define this as the direction of the chord in the last few microns as the 
track approaches the’ surface of the emulsion, the average error can be 
evaluated : 

From the formule given by Rossi and Greisen (1941), and using their 
notation we easily obtain for the angle between chord and trajectory 


(2) =4/3(t/w?), 


TABLE. IT, 
Values of pg/pg=5 55/g ; H=29,200 ; 1=1 cm. 


Mass of Particle 


Residual 

Tee 200 m, 300 m, 2000 m, 
be (per cent) | (per cent) | (per cent) 
100 2-7 2-9 4-3 
200 2:3 2-55 4-0 
400 2-0 2-15 3-4 
600 1 136 2-03 3-0 
800 LS 1-9 29 

1000 1-0 1-8 2-7. 


The calculations have been made for a particle of charge ¢, and masses. 
~200m,, 300m, and 2000m, for a path length in the air-gap of 1 cm. 
The values of py/ps, when multiplied by 2/(2—m) (about 1-4 in the present 
instance) give the error in the mass due to the scattering in the air. 


where t is the length of the chord in radiation units, and w is a function 

of the momentum and velocity. On the other hand, for the total angular 

deviation @ of a track from its original direction, after a path, t, we have 
(022 2t epee bi5 CE"). 

It follows that ¢€?) is easily evaluated from experimental values of (6°) 

(Goldschmidt et al. 1948). . 

A simple and more direct method of estimating the importance of 
all the above sources of error is to make an exposure without a magnetic 
field. The measured ‘“ deflections’? are then due to scattering or 
distortion. An examination of a plate containing a total of 50 tracks 
(10 mesons and 40 protons) has shown that the errors in the angular 
measurements are distributed in a gaussian curve, the half-width of 
which is of the order of 1-4°. 
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(iv) Spurious coincidences. 

The number of “ spurious coincidences ’’—the number of cases in which 
two unrelated tracks are attributed to the same particle—depends on 
the criterion adopted for the recognition of true pairs—and the method 
employed has been described in §2. Consider a track T, in P,, 
and let its angle of “‘ dip ” be ¢,+ 4¢, its residual range, 7,, and its direction 
x,. The region on P, in which a search for the partner track will be 
made, is the area R,, shown in fig. 3; this area is limited by two circles 
of radius 5 tan (¢,—4¢) and 6 tan (¢,+4¢) on the one hand; and by 
the circular paths C, and C,, corresponding to the maximum deflection 
expected for the particles with mass greater than 100 m.e., on the other. 

Suppose now that the tracks in P, are isotropically distributed. This 
assumption is only approximately correct, but it can hardly lead to serious 
errors. Suppose, further, that the tracks occur with a frequency per 
unit area equal to y,. The number of * fortuitous coincidences ”’ will then 
be given by the number of tracks contained in Rg, crossing the surface 
of the emulsion, with the correct angle of dip, 6+ 4¢, and correct direction, 
%o+4a. If dd is small, the required value is given to a sufficient 
approximation by 


2v Aa(dd)?.d3 tan? ¢ 
1= —, ——_“—_. : 
7 Pmin COs d 
For N tracks, the total number of “‘ coincidences ”’, C, is, 


N ‘ 
C=2n- 


1 
If the N tracks are isotropically distributed within the angles ¢, and dy, 
we then obtain 


co 2Nv Ac (4$)259 


Ps 
.| + tan? d6—tan +4] as 
TPmin E ? ? te 
In the conditions of our experiment : 
Au=1°=0-017 radians, y 
Ad=2°=0-034 radians, 
5=0°3 cm., 
Pmin =2 em, 
Considering only particles having a path, J, in air between 1 and 4 cm., 
we find, from Table III., that C~10-4vN. 


TaBe IIT. 
Values of @=[4 tan? d—tan +4] 


1 (em.) 1-0 1:5 2-0 3-0 4-0 
¢ (rad) ]-27 1:37 1-42 1-47 1-50 
@ (¢) 10-4 35:2 92:7 324:5 788-0 


If v< 100 cm. the number of fortuitous coincidences will be less than 
1 per cent. 
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In this calculation, it has been assumed that the correlation of the 
events has been made only on the basis of geometrical considerations 
whereas in actual fact, the observation of a close similarity in the prain- 
density of the two tracks of a pair has always been considered an Srential 
additional criterion. We have assumed, however, an isotropic distribution 
of the tracks which (see fig. 7) does not, in fact, exist. Approximate 
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Logarithm of radius of curvature p (cm.) versus logarithm of residual range r. 
The lines represent the function with slight corrections to allow for devia- 
tions of the range-energy relation from the power law for short ranges. 


estimates suggest that the two effects compensate one another, and it is 


reasonable to suppose that the calculations give a satisfactory indication 
of the order of magnitude of the number of spurious coincidences. 


§5. RESULTS. 
In fig. 4 the results of measurements on the tracks of 150 particles, 
observed in Kodak NT4 emulsions, have been plotted. The probable 
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error in the measurement of the radius of curvature, calculated in each 
case as discussed in §4 is indicated for each point. For the sake of 
clarity, the error in range, which is estimated to be smaller than 5 per cent 
in every case, has not been shown. A similar graph has been constructed 
for the other 210 measurements, made with Ilford G5 plates, for which 
the range-energy relation is slightly different, and the results are essentially 
similar. It will be seen that the results clearly show the division of the 
particles into groups corresponding to mesons, protons and heavier 


Fig. 5. 
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Mass distribution of 360 recorded particles. The dotted lines correspond to. 
mass values of u-mesons (215 m,), 7-mesons (285 m,), protons, deuterons 
tritons. It follows from equation (4) that a particle of charge Ze appears. 
(using this method) to have a mass of m/Z*, a being 0-59 for Kodak NT4 
and 0-56 for Ilford G5. a+ and o_ mesons are indicated in black. ; 


particles, the most prominent group being that due to protons. The linear: 
distribution, anticipated from equation (4), for particles of a given mass. 
and of range greater than 200, is clearly displayed. 

The results of measurements of the mass of 359 particles are shown 
graphically in fig. 5, and details of the measurements on the group of” 
particles with masses less than 400 m_,, including the secondary phenomena 
observed at the end of their range, are given in Table IV. 
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TABLE IV, 


Analysis of Mesons—Mass Measurements. 


9 


Plates Sensitive to Minimum Ionization Tracks. 


The mesons are classified phenomenologically, according to the secondary 


effects observed at the end of their range, into z, o, p (Powell, 1949) 


y-meson.is ene which is observed to originate by the decay of a z-particle 


in the emulsion. 


ow Residual Frobable : 
No Sign |, p Se Mass | errorin | Observations 
Deg.)|(em.)} Range Maas 
KO, ~12°)) — *| 13-0 | 1-94 105 220 20 p electron decay 
eee bate te 17-2} 2-18 108 240 30 p electron decay 
eal en AG me ws 5) 415 178 20 p 
-29 | + -| 29-5 | 2:94 308 215 10 p 
26 =. | 23°0+)'3-63 921 215 18 p 
76) -+- | 19-0:| 2-24 130 235 15 p electron decay 
DOT ew ol 182-80) 178 285 35 amp decay 
110 | 241 2-22 137 230 15 p (slow electron) 
158 30-2 | 2-71 344. 210 18 p electron decay 
164 |) + 37-5 | 2-56 258 220 10 be 
72 — | 42-0 |} 2-96 126 350 50 o (5 branch star) 
55:>| *— | 26-7) 2-76 | 327 216 10 p (blob)? 
125 61-5 | 3-63 650 215 i! p (slow electron) 
142 + {51-5 | 3-1 395 206 10 » (electron decay) 
39 | — |35-5| 3-84] | 249 390 30 | p (blob) 
K.2 24 — | 22-5 | 3-09 295 270 10 o (3 branch star) 
a 89} = (25-0 f2-54| 585 | 158 | 25 |p 
75 + | 21:0 | 2-3 57 240 30 7 
7 2 = 15:8 | 3-55 388 280 20 7 
6 “3 8-0 | 3-7 286 330 30 7 
96: (0 =- 433-0 12-5 305 190 30 p electron decay 
DS ee oO R27 315 120 18 p electron decay 
JS te 0 ES Da BS) 220 155 20 p electron decay 
141 — | 16:0 13-94} 690 250, 20 o (1 branch star) 
agp] eet 0-3'1'8-8 4 |" 933. °|.. 320 30 | p slow electron 
SOS oe 29'0' | 2-7 310 212 20 pw electron decay 
10| -+ | 10-0) 2-85} 385 | 220 20 | pw electron decay 
geile = 23°71 1-95) 145° | 205. | 20 7 
54 | — |444)1-85] 101 220 7 p (blob)? 
121 | ‘— | 10-:0|3-33| 357 265 30 p (blob) 
134); + | 12:0] 2:5 124 310 30 7 
190 | -—- | 13-0] 2-35) 235 200 40 p (blob)? 
191 | + | 11-0 | 2:7 206 250 20 
173 | + }12:0)2-:94] 335 225 15 p electron decay 
176.) 8:0 |} 2-88} 290 235 20 p electron decay 
177 | ‘= | 19-6 12:88 | 358 220 12 p electron decay 
1.8 64] — | 31-6 | 2-9 298 235 10 p slow electron 
Ga Te 1 aailio-98 || 165) 1320 W228. TT 
81 | + | 13-97 2-58] 270 210 25 pu electron decay 
85 — |13:5|3:55| 430 265 25 o (2 br. star) 
93 | + | 2441444] 1530 190 15 p electron decay 
91 + ‘+ 10-0 | 26 241 230 20 p electron decay 
99 = | 125 eo 200 » 165 25 p 
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TaBLE LV. (Cont.) 
peat "hig Wg an A ya 


a p |Residual Pavan Ob re 
v Si ss | error in servations 
No Sign Deg.)|(cm.)} Range Mass eee Z 
154 | — |22-5|2-78| 396 195 10 p (blob? and slow 
electron) 
164 | + |14-:0}2-14] 160 170 25 p electron decay 
174 | + | 16-5 | 3-0 370 225 15 p 
171 | + |14-0]2-05] 510 250 Bye ee 
140 | — | 19-0) 1-7 120 180 50 p 
226 12:0] 2:55] 140 292 50 o (1 heavy trk.) 
1.9. 29 | — |.17-0 | 2-75 202 260 20 p 
46} + |15:0]445] 714 280 20 7 
47 | + |17-5| 2-6 200 245 20 pe (incomplete) 
52 | + | 14-4] 3-0 400 215 5 p electron decay 
65 | + | 25:5] 2-45] 295 192 10 p electron decay 
74) — 9-0 | 4-0 450 300 20 p 
79 | + | 10-5] 2-35] 236 200 40 7 
95} — | 17-1) 4-0 508 290 20 p slow electron 
105}. — 9-0 | 3-0 108 350 40 a (1 heavy br.) 
121 | — }19-0| 3-33] 392 245 20 p (blob)? 
125 = 19-8 | 4-0 484 295 30 p (blob) 
3 | — |14:5]3-28] 496 218 25 p (blob)? 
Bol ate) 20:04" b-Oo' | © aes 207 15 p electron decay 
16 |. + 2-O1i3-7 460 260 20 ‘| -p electron decay 
30 | -— | 13-5). 4:5 902 250 40 p electron decay 
95 | + 9-0 | 4:75 | 1700 200 20 p electron decay 
84 = 9-0 | 3-2 227 310 30 o (4 branch star) 


Plates not Sensitive to Minimum Ionization Tracks. 


K.4. 26| — | 14:0] 2-34] 120 262 25 o (4 branch star) 
63 = Heer ip lioris: 160 160 30 p 
71 | = "4.19-0.1 2-99") 2 250 195 20 | p 

30 | — | 24:5] 2-99] 512 207 15 p (slow electron) 
31 — | 16-5} 2-81 267 240 25 p 
38 | + | 91-6 | 2-5 228 420 10 p 
41 — | 22-5 | 2-3 150 230 15 p 
44 + | 26:0 | 3-65 800 225 15 p 
91 — |13:-4]3-48] 800 208 15 p 
28 | — | 21-0] 3:48] 374 282 15 p 

107 — | 10-8.) 2-66), 220 245 30 o (3 branch star) 
47 |. --. |-20-0,).2:927) 9 346 230 20 p 
11] + 9-5 | 5:7 1300 340 80 7 
100 | + | 14:9] 4-72] 1030 292 20 1 
K.6, 68 : 21:5 | 3-26 205 330 30 p 
70 + | 17-0 | 2:88 197 283 7 7 
114 11-0}2-05] 85 250 50 | p 
92 eel: Ol eZ 45 258 20 7 
97 + |330| 2-17 91 262 13 7 
19 L 17-0 | 2:5 160 260 15 7 
26/ — | 285/31 | 450 228 nam p 
45 — 30-0 | 2:07 153 200 17 p 
53} + |12-5)3-45| 593 235 16: |p 
OT _ 11-0} 2-94 160 315 40) p 
109 | — | 26-4 | 2-5 239 217 10 p 
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The observed distributions of the values of the mass of the particles, 
taken separately according to their sign, are shown in figs. 5(b) and (c). 

The distribution in mass of the positively charged particles (fig. 5(b)) 
shows : 

(a) A peak due to particles with masses in the interval from 160 to 
250 m,, the mean value being 217-4 m,*. About 80 per cent of these 
particles decay with the emission of an electron, of which the track can 
be distinguished. Six of them are y-mesons created by the decay of 
m-particles in the emulsion. 

(b) A peak, just resolved from the first, and due to particles with 
masses between 240 and 340 m,, the mean value being 281-10 m,. All 
these particles decay at the end of their range with the emission of a 


p-particle. 
(c) A more prominent group, corresponding to the particles with a 
mean mass of 1870m,, +27m,. None of these particles give any 


indication of inter-acting with matter at the end of their range, nor do 
they produce any observed charged particles by spontaneous decay. 

(d) A peak of lower intensity can be seen at about 3700 m,, and an 
indication of a broad distribution due to particles of a greater mass. 
None of these particles give any indication of inter-acting with matter, 
nor of undergoing spontaneous decay. 

The distribution of values of the mass of negative particles (fig. 5(c)) 
shows : 

(ec) A peak between 160 and 250m,, and a broader distribution 
extending up to 400 m,. Of the particles with masses between 160 and 
250m,, only four show electron decay; and five others give rise 
to slow electrons originating at the end of their range. Of the 
particles between 250 and 400m,, 9 produce nuclear disintegrations, 
and the others stop without leading to the production of any secondary 
charged particles for which the evidence is unambiguous. Jn ten cases 
there is an unusual thickening at the end of the track, see PI. II.(5), 
due to an assembly-of two or three grains which might be attributed 
to a secondary particle of very short range ; but in all cases except three 
it has been impossible to decide whether this ‘‘ blob” was produced by a 
slow electron or by a heavy particle. Inthe three cases the “ blob ”’ was 
definitely too thick to have been produced: by an electron. 

(f) In fifteen cases, a particle which crosses the air-gap and is brought 
to rest in one emulsion, originates in the emulsion of the other plate. 
The complete track therefore has two “ends”. The direction of motion 
of the particle, in events of this type, can be determined by the increased 
grain-density in the track as the particle slows down, and the more 
numerous deviations due to Coulomb scattering towards the end of the 
range. The particles of this type are all positively charged and the 
distribution in the observed values of their mass is shown in fig. 5 (a). 


* The errors, o, are calculated from the residuals, v;, according to the formula 


c 0? 
°~ A} n(n—I)* 
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(g) The angular distribution, relative to the vertical, of the directions 
of motion of all particles on which measurements have been made, 
at their points of entry into the first emulsion through which they pass, 
are shown in fig. 7. In this figure, the particles have been grouped 
together according to their masses and the signs of their charge. In 


Fig. 6. 


negative particles. 
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Effects of distortion of the gelatine on the precision of the measurements. 
(a) 61 tracks from 2 pairs of 50» thick emulsions (25 days exposure). 
(b) 56 tracks from 1 pair of 100 thick emulsions (15 days expos-re). 


(c) 69 tracks from 1 pair of 200 thick emulsions (10 days exposure) 
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(b) w-; (¢c) a-; (ad) at; (e) recoil particles ; 


(f) protons, deuterons and heavier particles. This histogram also contains 
some tracks which were not in a suitable position for being measured, 
although their identity could be established with certainty. 
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most cases there is a tendency for the directions of motion to be grouped 
round the vertical, but to be deflected from it in a sense depending on the 
sign of the charge. ; 

(h) There is a gradual reduction in precision as the thickness of the 
emulsion employed is increased—see fig. 6. 


§6. INTERPRETATION OF RESULTS. 

From the evidence given in sections (a) to (g) above it is concluded that : 

(«) The positive particles corresponding to group (a) are j.+-mesons 
created by the decay of z+-particles. In those cases in which the 
u*-particle has an energy greater than 4:2 MeV., we must assume that the 
parent 7'-particle decayed whilst in motion. The average value of the 
mass of these particles is 217+4m,. We therefore identify them with 
the mesons of the penetrating component of the cosmic radiation, of which 
the mass, as determined at Berkeley, is 215+4m,. 80 per cent of 
these particles decay with the emission of an electron. In view of the 
fact that “thin” tracks inclined steeply to the plane of the emulsion are 
difficult to detect, and since there is an isotropic distribution of the 
directions of motion of the decay electrons, we may assume that we fail 
to detect 20 per cent of the electron tracks, of which the specific ionization 
has the minimum value, in the conditions of the experiment ; but that, 
in fact, all the w*-mesons decay with the emission of an electron. 

(6) The particles described in group (6) are obviously 7+-particles. 
From the observed value of the mass, 281+7m,, they can be identified 
with those produced artificially. 

(y) The negative particles described in group (e) are made up of 
p- and z--particles. Nine of them produce observed disintegrations, 
and of these the average value of the mass is 288-++13m,. They can therefore 
be identified as z~-particles. The average mass of the other 38 negative 
particles, only four of which appear to decay with the emission of an 
electron, is 236+8m,. It is reasonable to attribute this value, inter- 
mediate between the masses of the 7 and yw particles, to the presence 
among the ~-mesons, of some 7~-particles which produce no visible dis- 
integrations. The observed value of the average mass suggests that 
28-44 of these mesons are .-particles, and 10-/4 are 7 -particles. 

(8) The group with apparent masses between 1100 and 2400 m,., which 
contains 48 per cent of all tracks measured, is formed by positive particles 
which stop in the emulsion without showing any inter-action with matter, 
nor any evidence of spontaneous decay. Their average mass is 
1870--27 m,, and they may be identified with confidence as protons. 

(c) The group between 2400 and 4400 m_,, of smaller intensity than. the 
protons and less clearly resolved, is also made up of positive particles 
which produce no secondary effects; they are therefore regarded as 
deuterons. | 

No accurate discrimination of particles of different types is possible 
in the region of the spectrum beyond 5500 m,._ This is due not only to 
the small deviations suffered by heavy particles, but also to the fact that 
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their charge, which may be different from |e], is unknown. We therefore 
attribute this group to a mixture of tritons and «-particles, which cannot 
be resolved in the conditions of our experiment. 

(C) Of the seven particles which appear to have an intermediate 
value of the mass, six are positively charged, and one appears to be 
negatively charged. We must expect that a certain proportion of the 
different particles which we record will be scattered through an 
angle of a few degrees or more, in passing through the gas between the 
plates, as the result of elastic collisions with oxygen or nitrogen nuclei. 
Such a deviation in the trajectory will usually lead us to reject the particle 
from among those which give tracks acceptable for measurement, for 
the directions of motion—at the points of exit from and entry into the 
two emulsions—will not in general lie on a circle. In a small proportion 
of such cases however, the deviation will be of such a character that the 
tracks are still acceptable. This will be so, for example, if the nuclear 
collision occurs near the mid-point of the trajectory of the particle in the 
air-gap. The observed deviation will then be less or greater than that due to 
the magnetic deflection alone, and the calculated values of the mass 
will be correspondingly greater or less than the true value. 

Other processes, such as inelastic scattering of protons in,the emulsion, 
will tend to lead to an over-estimate of the masses of the particles. In 
addition, as has been shown in §4, we must expect to observe a certain 
number of fortuitous coincidences, and their number in the particular 
conditions of the experiment should be 3 or 4. Because of these and 
other possible processes, it appears reasonable to attribute the observation 
of particles.with intermediate values of the mass to scattering in the gas, 
or in the surface layers of the emulsion, or to other uncontrolled sources 
of error. The results are therefore regarded as giving no support for the 
view that other particles exist which are sufficiently stable to be brought 
to rest in solid materials, and which have masses in the interval between 
7-mesonsand protons ; and as proving that if such particles exist they either 
occur in matter bombarded by cosmic radiation at 11,000 feet, with a 
frequency less than 3 per cent of that of 7 and yw particles and less than 
1 per cent of protons; or, alternatively, they have a lifetime shorter 
than 10~-?° sec. 

(yn) Of the particles which originate in the emulsion (see fig. 5(a)), 
11 have masses equal, within the limits of error, to that of protons. 
Three of the others may be identified as deuterons and one is of greater 
mass—probably an «-particle. The majority of the particles may be 
attributed to knock-on protons produced by the impact of neutrons 
with hydrogen nuclei ; and others to nuclear disintegrations produced 
by neutrons of moderate energy, in which a single charged particle is 
emitted. 

(9) The distribution in the directions of motion of the pt and ~~-mesons 
—see fig. 7 (a) and (b)—shows that most of these particles move down- 
wards at a small inclination to the vertical —see (Camerini ct al. 1948), 
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There is some indication that the maxima in the two distributions 
deviate from the vertical in opposite senses according to the sign of the 
charge. This may be accounted for if a considerable proportion of the 
particles are generated in the atmosphere or other material above the 
apparatus. They will then travel in the magnetic field for a distance of 
several centimetres before entering the emulsion, and will suffer a 
corresponding deflection. 

The statistical weight of the angular distributions of the 7+ and 
a7 -particles is too small for any conclusions to be drawn, but the 
observations are consistent with the assumption that the majority of 
particles are locally generated in the matter in the immediate neighbour: 
hood of the plates, and that the directions of motion are distributed 
approximately isotropically. 


TABLE V. 


Type of | Number of Number with Percentage of 


yarticle | particles all particles 
: i electron decay | slow electron “blob ’ 


ut 30 19(24)* =m Aa 8541-5 
Ve 9844 4(18)* 6(18)* 6(18)*2 84125 
7 14 = Ba ee 4-0 +1-0 
ae 19+4 Jt peri) 3 5-4 42-0 
H} 168 zt ses = 4844 
H2 55 = = = 1642 
H3,a etc} 39 ee be Be 1142 


* These values are deduced from observations on electron-sensitive emulsions 
only. The figure in brackets gives the corresponding total numbers of particles 
observed in these plates. 


The distributions for the protons, fig. 7(e), can be regarded as being 
made of two groups, one of which is due to particles moving in directions 
inclined to the vertical at less than 40°. From the displacement of the 
of the peak from the vertical, 10°, we must assume that most of the particles 
have travelled at least 5cm. in the magnetic field. The other group, 
containing about 60 per cent of the total, is made up of particles of which 
the directions of motion are more widely dispersed, and they may be 
attributed to protons which originate in nuclear explosions produced 
in the matter of the apparatus. 

The results of the present experiments are summarized in Table V., 
which shows details of the proportions of the particles of different types 
amongst all those on which measurements have been made. 


§6. COMPARISON WITH OTHER METHODS. 


A comparison of the present results on the masses of the particles 
with those obtained in other experiments employing photographic 
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emulsions is shown in fig. 8. The figure gives the mass distributions 
obtained by Goldschmidt et al. from scattering measurements, and by 
Lattes, Occhialini and Powell from “grain-counts”. In each case ‘it 
can be seen that the “width” of the distribution, due to particles of a 
single type, is much greater than that obtained by the magnetic deflection 
experiments. 
Experiments similar in principle to those described in the present 
paper have been carried out by Alichanian and Alichanov and their 
collaborators (1948). Using Geiger—Muller counters, they have measured 
the deflection, produced by a magnetic field, of particles brought to rest 
in a certain thickness of lead. The mass distribution thus obtained 
shows, (a) a prominent peak corresponding to the protons ; (6) groups of 
particles of greater mass than protons, but not individually resolved ; 


Fig. 8. 
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Comparison with other methods. (a) The mass distribution of all the p-mesons, 
(7. e. w+, w~ and that fraction of 7~ which does not produce stars) obtained 
by magnetic deflection (dotted line), and that obtained by scattering 
measurements (full line) (Goldschmidt ef al. 1948). (b) The mass 
distribution of the ~+ obtained by magnetic deflection (dotted line), and 
that obtained by grain-counting (Lattes et al. 1948). 


and (c) a continuous distribution towards lighter masses, which fills 
completely the interval between the groups due to protons and u-mesons. 
The Soviet experimenters have interpreted their results as indicating 
the existence of many types of particles with intermediate masses which 
can carry both positive and negative charges. Further, a preliminary 
report has been given by Brode (1949) of experiments of the same type 
but capable of greater precision, in which the deflections are measured A 
observations with Wilson chambers. In this case, it has been suggested 
tentatively, that the observations may correspond to the existence et 
positive particles of mass about 1000m,. In the case of both these 
experiments the particles are arrested in layers of lead, and it is therefore 
not possible, in general, to observe any secondary processes occurring 
at the end of the range of the particles. 

The conclusions to which we have been led by the results of the present 
experiments are therefore in contradiction with those obtained by other 
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methods. The seven observed particles with intermediate values of the 
mass have been attributed to {incontrolled sources or error. It appears 
to be of great significance that all but one of the particles are positively 
charged, and that they show no evidence of spontaneous decay at the end 
of theirrange. A similar feature is characteristic of the results described by 
Brode, inwhich all the observed particles of intermediate mass are positively 
charged. In view of the great variety of physical processes which can 
occur in apparatus of the type described by Alichanian and his colleagues, 
in addition to those which alone are assumed to be operative in the 
interpretation of the results, it may be suggested that the existence 
of other types of mesons, in addition to the 7- and j-particles, should be 
accepted with very great reserve; and that the present experiments 
prove that if they exist, they are certainly very much more rare than 
. previous experiments have been taken to indicate. 
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Puates I. & IL. (a). 
Examples of “ partners ”’. 


Photographs of the tracks of particles at their exit from the first emulsion 
(first row), at the entry into the second emulsion (second row), and at the 
end of their range. ; 


The great similarity between “* partners ’’, which is immediately apparent 
under the microscope, has been partially lost in the photograph. 


Puate IT. (5). 


Negative meson. The arrow points to the “ Blob.” at the end of the range, 
which seems to be characteristic of negative mesons only. 
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VIE. Notices of New Books and Periodicals received. 


Supersonic Flow Around Cones of Large Yaw. (Massachusetts Institute of 
Technology Technical Report No. 5.) [Pp. 125.] (Cambridge : Massachu- 
setts. ) 


Tur properties of the flow about a yawing cone can be obtained as a power 
series in the yaw ¢ ; earlier reports in this set (Nos. | and 3) presented tabular 
results for the coefficients of the terms independent of ¢ and linearly dependent 
one. The present report (in which the same value 1-405 for y is used) deals 
with the ccefficients of <? for cones of semi-angle between 5° and 25° 
[5(2-5)15, 20, 25] and gives in Part I tables of the additional second order eftects 
on. velocity, pressure and density. Part II contains information for deriving 
drag (which has a second order contribution in e), normal force, ratio of shock 
wave to cone yaw, and the ccefficients determining the variations in, entropy 
and in the shape of the shock wave which to this order ceases to be strictly 
circular in cross section. [new 


The Origins of Modern Science. By Herpert BurrerFreLp. (Bell & Sons.) 
Price 10s. 6d. net. 


Tas work by the Professor of Modern History in the University of Cambridge 
is a welcome contribution to scientific literature. 

A number of eminent scientists by their writings have demonstrated the 
cultural aspects of scientific thought and experiment. It is now gratifying to 
record the entry of a professional historian into this field in an endeavour to 
build a bridge between the humanities and the sciences. 

In the reviewer’s opinion the writer has achieved a conspicuous success. 
For this reason this is not the place to make minor criticisms. When the 


average student of science is so well versed in the history of his subject that he 


is capable of crossing swords with Professor Butterfield on his interpretation of 
some point in the views of, let us say, a sixteenth century physicist, then the 
bridge that the author is seeking to build will be truly fashioned. For this 
reason alone this volume is particularly to be commended to those who in 
Universities are concerned to see some broadening of the training of the 
specialized scientist. 


An Advanced Treatise on Physical Chemistry.—Vol. I. Fundamental Principles 
and Properties of Gases. By J. R. Partineron. [Pp.943+-40.] (Longmans. ) 
Price 80s. 


Turis book, as the author points out, is advanced in the sense of being 
comprehensive, and not in the sense of assuming much previous knowledge ; 
for example, no previous knowledge of trigonometry or calculus is assumed, a 
concise account of these subjects being given in the mathematical introduction. 
This introduction is followed by sections on thermodynamics, kinetic theory of 
gases, statistical mechanics, quantum mechanics, thermometry, and a long 
final section on the properties of gases. Most sections start with a 
brief historical introduction and include a detailed account of the contributions 
of various workers in the field, with tables of experimental data and numerous 
references; the latter number about 18,000. Emphasis is laid throughout on 
the experimental side, though the theory is also thoroughly described. This 
book will be found a very valuable work of reference. tea aa AE 
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Applied Mathematics for Engineers and Scientists. By S. A. SCHELKUNOFF. 
[Pp. 472.] (D. van Nostrand Company Inc., 1948.) Price 36s. 

PRESENT day science demands a mathematical discipline which often cannot be 
obtained from an undergraduate course. The author aims at providing this 
background for students of applied science. An idea of the ground covered 
gan be obtained from the following list of contents: Complex numbers— 
Theory of approximation—Solution of equations—Power series—Difierentiation 
—Integration—Vector analysis—Coordinate systems—Exponential functions— _ 
Differential equations of the first and higher orders—Partial differential 
equations—Conformal transformations—Contour integration—Linear analysis 
—The Gamma function—Exponential and Fresnel integrals—Bessel and 
Legendre fanctions—Formulation of equations. Thus the first half treats the 
general mathematical methods and the second half introduces the frequently 
occurring transcendental functions. The portion dealing with differential 
equations is particularly well written. The introductions to the use of Green’s 
function and the Laplace transform method are two of the most valuable 
features of the book. The principles involved have been illustrated by concrete 
examples as far as possible. A chapter on tensors and matrices would perhaps 
have increased the general usefulness of the book. On the whole the book 
deserves to be widely used. da hy, bea es 


The Theory of Atomic Collisions. By N. F. Morr and H. 8. W. Massey. 
Second Edition, Royal 8vo. {Pp. 388.] (Oxford: The University Press, 
1949.) Price 35s. net. 


In this work, which has been a classic since the publication of the first edition 
in 1933, the authors set themselves the task of discussing those collisions which 
do not involve creation or annihilation of particles, and can therefore be dealt 
with by wave-mechanics alone. In the second edition the form of the book 
remains unaltered, but the subject-matter is extended to include the phenomena 
of classical nuclear physics. There are valuable references to a large number 
of papers on, both atomic and nuclear problems published since 1933, and the 
original treatment is expanded at varous points. 

Broadly speaking, the book is divided into three parts. The first eight 
chapters are devoted to the general theory of collision processes, and include 
subsidiary mathematical and physical development where necessary. The 
main, addition here is in the treatment of the dispersion formula, and of the 
collision complex or compound nucleus, although a more detailed analysis of 
the Born approximation has also been made. 

Chapters IX. to XIII. deal with applications to atomic and nuclear collisions. 
More stress is now laid on nuclear problems and the chapter on the collision of 
electrons with molecules has been omitted. 

In the last chapter relativistic two-body problems are considered. In 
orthodox quantum mechanics it has been usual to treat such problems by 
quantising the radiation field, thus bringing them outside the scope of this 
book, but the alternative method of the retarded potential has also been 
available, and the authors use this instead. Recently Feynman has criticised 
the use of field theory and developed the second method further, but it was 
obviously impossible to include a discussion of his work here. Some important 
work by Stueckelberg (1938), who foreshadowed Feynman’s approach, has, 
however, apparently been overlooked by the authors. This chapter now includes 
a brief treatment of certain meson problems, together with useful references. 

Altogether this new edition should prove invaluable to a wide class of 
physicists. Kiv. Re 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 


